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Abstract: We investigate the large N instanton effects of partition functions in a class 
of AA = 4 circular quiver Chern-Simons theories on a three-sphere. Our analysis is based on 
the supersymmetry localization and the Fermi-gas formalism. The resulting matrix model 
can be regarded as a two-parameter deformation of the AB JM matrix model, and has richer 
non-perturbative structures. Based on a systematic semi-classical analysis, we find analytic 
expressions of membrane instanton corrections. We also exactly compute the partition 
function for various cases and find some exact forms of worldsheet instanton corrections, 
which appear as quantum mechanical non-perturbative corrections in the Fermi-gas system. 
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1 Introduction 

Over the last couple of years, many interesting features of non-perturbative effects in Id- 
theory on AdSii background were discovered via the AdS/CFT correspondence. Utilizing 
the powerful techniques of the supersymmetry localization [1] and so-called Fermi gas 
approach [2] to the partition function of the dual CFT on 5^, now we have a very detailed 
understanding [3] of the non-perturbative effects in M-theory on AdS^ x S'^ j'Lk-, which is 
holographically dual to the ABJ(M) theory [4, 5]. 

It is realized that the existence of the two types of instanton effects, worldsheet instan- 
tons [6] and membrane (D2-brane) instantons [7], is necessary for the consistency of the 
theory. From the bulk M-theory perspective, these two types of instantons are both origi¬ 
nating from certain configurations of M2-branes wrapping some three-cycles, but they are 
distinguished by the different dependences on the Planck constant h = 2TTk of the Fermi gas 
system. Membrane instantons are already visible in the semi-classical small h-expansion 
of the Fermi gas, while worldsheet instanton effects are non-perturbative in h. One can 
study such worldsheet instanton effects from the opposite large h regime, by 1/h-expansion 
of the matrix model associated with the Fermi gas system, which corresponds to the or¬ 
dinary genus expansion of the matrix model with the string coupling given by = 1/h 
[8-10]. From the viewpoint of this g^-expansion, the membrane instantons appear as non- 
perturbative effects in gg- The pole cancellation mechanism found in [11] is important for 
the consistency of these two expansions, since it guarantees that we can go smoothly from 
the small h regime to the large h regime. There are also additional contributions, namely 
the bound states of worldsheet instantons and membrane instantons [12], which are hard to 
study from both small h and large h expansions. Fortunately, in the case of ABJ(M) theory, 
we have a complete understanding of the worldsheet instantons, membrane instantons, and 
their bound states, thanks to the relation to the (refined) topological string on local x P^ 
[3, 13-16] and exact results for various specific values of the parameters [13, 14, 17, 18]. 
(see [19, 20] for similar progress in half-BPS Wilson loop) 

However, for more general 3d gauge theories with less supersymmetry, we still do 
not know detailed structures of the non-perturbative effects^. Some progress along this 
direction has been made in the study of an AA = 4 7/ {N) gauge theory with Nf fundamental 
and one adjoint hypermultiplets, which appears as the worldvolume theory on N D2- 
branes in the presence of Nf D6-branes. After applying the localization technique [22- 
25], the partition function of this theory on is reduced to a matrix model, called the 
Nf matrix model [26, 27]. In [28], using the Fermi gas approach with the identification 

^The only exception so far is the orbifold ABJM theory analyzed in [21]. The grand potential of this 
theory has a simple relation to the one of the ABJM theory. 
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h = Nf, the analytic forms of the first few membrane instanton corrections of this model 
were determined. Worldsheet instantons can also be studied, in principle, by the genus 
expansion of the Nf matrix model. The genus-zero and the genus-one free energies of the 
Nf matrix model were calculated in [26], but in practice the computation of the higher 
genus corrections is not so easy. Instead, the analytic forms of the first few worldsheet 
instantons were found in [28] from the exact computation of the partition functions at 
finite N up to certain high N = A^max- It turned out that the results of the membrane 
instantons and worldsheet instantons in the Nf matrix model are quite different from the 
ABJ(M) case. In particular, the membrane instanton coefficients are given by the gamma 
functions of Nf and quite different from the Gopakumar-Vafa type formula [29] in (refined) 
topological string, where only trigonometric functions oi h or 1/h appear. To understand 
the underlying structure better, it is desirable to study non-perturbative effects in various 
other models with AA = 4 supersymmetries. 

In this paper we study a special class of AA = 4 super Chern-Simons matter theo¬ 
ries in three dimensions [30, 31]: a circular quiver gauge theory with the gauge group 
U{N)k X U{N)q~^ X U{N)-k X and bi-fundamental hypermultiplets one by one 

between nearest neighboring pairs of gauge groups. The subscripts in the gauge group rep¬ 
resent the Chern-Simons level for each factor. In this paper, we will refer to this theory as 

9 )"aiodel”. It is expected that the (p, g)-model is the low-energy effective theory of N 
M2-branes probing the orbifold ji^p x Zg))/Zfc [32, 33], where the orbifolding action is 
given by (2.90), and this model is dual to M-theory on AdS^^ x {S'^l{'Lp x Zq))/Zfc through 
the AdS/CFT correspondence. The (p, g')-model can be regarded as a two-parameter de¬ 
formation of the ABJM theory, hence it is expected that there is a rich non-perturbative 
structure in this model. For instance, the ABJM model corresponds to the {p,q) = (1,1) 
model, while the Nf matrix model corresponds^ to the {jp,q) = (1,A^/) model with the 
Chern-Simons level k = 1. 

We will study the large N instanton effects in the (p, ( 7 )-model by analyzing the par¬ 
tition function on . By applying the localization method, the partition function of the 
(p, g)-model is reduced to a matrix integral [22-24], which can be further studied by the 
Fermi gas formalism with the identification h = 27rk. Note that the partition function of 
the {p, ( 7 )-model is invariant under the exchange of p and q with fixed k. In the original 
set-up of the circular quiver gauge theory, p, q and k are all integers, but at the level of 
the partition function, we can consider an analytic continuation of the parameters {p, q, k) 
to arbitrary real (or complex) numbers. The study of the {p, (;)-model in the Fermi gas 
formalism was initiated in [37, 38]. The perturbative part of the grand potential was de¬ 
termined in [37], and in [38] it was found that there are three types of membrane instanton 
corrections in the grand potential, 

Oie-N), 0(e-^), Oie-f^), (1.1) 

where p denotes the chemical potential of the Fermi gas system. These instantons con- 
tribute to the canonical partition function at large A^ by 0(e V p ), 0(e V <i ) and 
^This correspondence is understood from the 3d mirror symmetry [34, 35] as explained in [4, 36]. 
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_ / pqkN 

0{e 2 respectively. The first two types of instantons are simply related by the 

exchange of p and q. In this paper, to study such membrane instanton corrections we will 
develop a systematic method for the h-expansion (WKB expansion) of the grand potential. 
Using the data of the WKB expansion, we determine the analytic form of the leading men- 
brane instanton correction of the first two types in (1.1) for generic (p, q, k), and for the last 
type in (1.1) we find the analytic forms of the leading and the next-to-leading menbrane 
instanton corrections. Also, there are worldsheet instanton corrections of the order 

( 1 . 2 ) 

-2tt 

which contributes to the canonical partition function by 0(e V ) at large N. From the 
exact computation of the partition functions at finite N, we find the analytic form of the 
leading worldsheet instanton correction for generic (p, q, k). We will see that our results of 
the menbrane instantons and worldsheet instantons satisfy the pole cancellation conditions 
as expected, and they also correctly reproduce the known results of the ABJM model and 
the Nf matrix model by taking the appropriate limits of the parameters (p, q, k). Study of 
the bound state corrections in the (p, q) model is beyond the scope of this paper. 

This paper is organized as follows. In section 2, after introducing the Fermi gas for¬ 
malism of the (p, ( 7 )-model, we explain our method of the WKB expansion of the grand 
potential and our algorithm for the exact computation of the partition functions. We also 
comment on the instanton effects in the (p, q) model from the dual gravity viewpoint. 
Then, we present the results for the (p, q) = (1, q) case in section 3, and the results for the 
general (p, q) case in section 4. In the both cases, we find the analytic forms of the first few 
membrane instanton and worldsheet instanton corrections. In section 5, we consider some 
interesting cases with the special values of (p, q). We discuss that the grand potential for 
{PiQ) = (1) “1) is captured by the refined topological string on a resolved conifold. We also 
find an exact closed form expression of the grand partition function of the (p, q) = (2, 2) 
model at A: = 1. Section 6 is devoted to conclusions and discussions. In appendices A to 
F, we summarize some useful results used in the main text. 


2 Fermi-gas approach to AA = 4 quiver CS matrix model 

In this section we introduce the ideal Fermi gas formalism of the (p, g)-model and discuss 
how to extract information on the large N non-perturbative effects in the corresponding 
M-theory dual. 


2.1 Prom matrix model to Fermi-gas 

It is known that the partition function of the (p, g)-model on is reduced to the matrix 
integral thanks to the supersymmetry localization [22-24], 


Z{N,k) = 


(iV!)P+'? J (27r) 


. P+9 

n 




p+q [)[ 


\N 




i<j 


2 sinh 


(a) (a) 2 

Pi -Pi ^ 


i=i 


n*j2cosh 


(a) (a + 1) ’ 

Ml -/^i 

( 2 . 1 ) 
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where . This matrix model can be further simplified by using the ideal 

Fermi gas approach [2]. To be self-contained, here we briefly review the derivation of the 
Fermi gas representation. First, by using the Cauchy determinant formula 


n 


i<j 


2 sinh 


2 sinh 


i V.i 


nij2cosh 2 

we rewrite the partition function as 


E(-irn 


j 2cosh^^^^’ 


1 


• p+i 


/ n 




N 


N\ 

(7 

p+q-l 


a=l ^ > j=l 


X 


n 


N 

n 


(a) (a+1) 

«=i b=i 2cosh^^^^^^- 

N 


i=i 2 cosh — 2 fc^ 


t7 j = l 

where we have trivialized most of the permutations and the function p is given by 


( 2 . 2 ) 


p{x,y) = 


2'Kk 




2 _ife / (q+l )\2 p+q-l 

e4:7vk^ 47r''^ ^ 


a =2 2 cosh 




a=2 


2 cosh 


fj_(a)_^(a+l) 

2fc 


2coshii5^' 


(2.3) 

Thus we can regard the partition function as an ideal Fermi gas system with the density 
matrix p(x,y). The expression of p is further simplified if we regard p as the matrix 
element of the quantum mechanical operator p{Q, P) with {Q, P) satisfying the canonical 
commutation relation 

[Q,P]=ih, h = 2'Kk. (2.4) 

Then, the density matrix p(x, y) is understood as^ 


p{x,y) = ^(x|/5|y). 


where 


KQ+) = 


— - 


2 cosh ^ 


-e 


2 cosh ■ 


p- 


(2.5) 

( 2 . 6 ) 


iQ^ 


Using the convenient formula e s'* f{P)& 2'* = f{P — Q) and e zr g{Q)e 2 h = g(^Q + p)^ 
we find 


_ i_ p2 

‘-tti X 


P{Q,P) = e 2 ft 


1 


1 


2 cosh $ I f^2cosh4 


i p2 

7e2ft^ . 


(2.7) 


®We are using the following convention 

{Q\P) = I ^\Q){Q\ = 1, I ^\P){P\ = 1. 
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If we perform the similarity transformation 


2 cosh 


O \ p2 ^_^ p2 / O N 9/2 


2/ 


pe 


2 cosh — 
2 


then the operator p is simplified to 


KQ,P) = 


1 


2 cosh ^ I 


qV/2 


2 cosh y 


2 cosh 


g/2’ 


( 2 . 8 ) 


(2.9) 


Note that the partition function is invariant under any similarity transformations. In the 
coordinate representation, the density operator (2.9) is expressed by 


pixi,X2) 


1 _ 1 _ 1 b(- + ~ ^ 2 ) P _ ijxi - 3 ^ 2 ) \ 

(2coshf)®/^ (2coshf)''/^ V2 27rA: ’2 2Trk J 


( 2 . 10 ) 

where B(x, y) = r(a:)r(y)/r(x + y) is the Euler beta function. For p S Z>o, one can show 
that this expression reduces to the one'^ in [39]. In particular, for p = q = 1, it reduces to the 
density matrix in the ABJM Fermi-gas as expected. Also, the case for {p,q,k) = (1,2,1) 
gives the density matrix for the U{N )2 x U{N +1)-2 ABJ theory [13, 40, 41] as explained 
in [28]. Note that the partition function is invariant^ under the exchange p q. This 
invariance is no longer manifest in the coordinate representation (2.10). 


2.2 Grand canonical formalism 

As was proposed in [2], a useful way to treat this system is to introduce the grand canonical 
partition function 

00 

E{K,k) = l+^K^Z{N,k), K = e^, (2.11) 

Af=l 

where k and p is the fugacity and the chemical potential, respectively.® We can return to 
the canonical partition function by 


The grand partition function is given by the Fredholm determinant 


( 2 . 12 ) 


E{k, k) = det(l + Kp) = ]^ (1 + ^An), (2-13) 

n=0 


^More explicitly, it is given by 


p{xi,X2) 


1 _1_1_1_ \ ( (2j-lp 

2{p-l)!,rfc (2cosh;^)‘*/2 2cosh21n^ ib=l 2^k J ^ 4 

1 _1_1_ XT-x,, TT^~^ \( x-,-x^ \2 . . 2 ] 

(2cosh;^)'>/2 (2cosh;^)'>/2 2sinh2U_^ = l 2-xk J \ 


for odd p 
for even p 


® We can easily show this by the canonical transformation (Q', P') = (— P, Q) and the similarity trans¬ 
formation (2cosh ■^)‘*/^(2cosh ^)P/^p(2cosh ■^)“'’/^(2cosh 
®In the following, we will use both re and p, interchangeably. 
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where (n £ ^>o) are the eigenvalues of p. The spectral problem for this system is 
represented by the Fredholm integral equation: 


/ CX) 

dx p{x,x')4>n(,x') = Xn(j)n{x). 

-CX) 

It is easy to see that the grand potential is given by 


OO /_ \l 

J (k, k) = log H(k, k) = 

e=i 


where Cp('S) is a spectral zeta function dehned by 


(2.14) 


(2.15) 


(2.16) 


(^{s) = Ttp^ = J2K- 

n=0 

In particular, for s = £ ^ ^>o, h can be computed directly by the multi-integral 

/ CX) ^OO 

■ dxi-■ ■ dxi p{xi,X 2 ) ■ ■ ■ pixi,xi). (2.17) 

-CXD J —OO 

As in [42, 43], it is convenient to rewrite (2.15) as the following Mellin-Barnes like integral: 

(2.18) 


rc+ioo 


= - / ^r(s)r(-s)Cp(s)e^^ 

J c—ioo 


The constant c must be taken such that 0 < c < 1. Depending on the sign of p, one 
can deform the integral contour in the following ways. For p < 0, one can deform the 
countour by adding an infinite semi-circle as shown in hgure 1. Then the integral can 
be evaluated by the sum of the residues over all the poles in the region Re s > c. As shown 
in [42], the spectral zeta function Cp(s) does not have any poles in the region Res > 0. If 
0 < c < 1, the poles inside the contour are located at s = .^ (.^ G ^>o)) coming solely from 
the factor r(— s), and thus we precisely recover the sum (2.15). On the other hand, for 
p > 0 one can deform the contour by adding the opposite semi-circle C- as in hgure 1. 
In this case, one needs the information of the poles in the region Res < c, in which Cp('S) 
may have non-trivial poles, and the problem is highly non-trivial for general k. As we will 
see below, in the semi-classical analysis, we can hnd all the poles of Cp('S) and compute the 
large p expansion systematically. 


2.3 Classical approximation 

Since the ideal Fermi gas formalism provides us with the quantum mechanical description 
of the system, it is useful to consider the semi-classical /i-expansion, or equivalently the 
small-fc expansion 

. OO 

JwKB{p,k) = (2.19) 

n=0 

where Jwkb is expected to approximate the exact grand potential J' up to exponentially 
suppressed corrections in h. Let us start with the classical approximation, namely 0{k~^). 
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-X- 



M -X-X^ 
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Figure 1. One can deform the integration contour in (2.18) into a closed path by adding an infinite 
semi-circle C+ {C-) for /i < 0 (/r > 0). We schematically show the poles of the integrand. 


Note that the large /r expansion in this limit has already been analyzed in [38]. Here, we 
simply re-derive their result by using the Mellin-Barnes integral (2.18). In the classical 
approximation, the density operator is given by 

“ (2coshf)«(2cosh|)"- 

Then, the spectral zeta function can be easily computed by the phase space integral 


(0) ^ dQdP 1 _ 1 ^ 1 r2(f)r2(f) 

^ i_oo 47r2 (2 cosh (2 cosh T{ps)r{qs) 

We would like to know the large expansion of the classical grand potential. 
(2.21) into (2.18), one obtains 


( 2 . 21 ) 

Plugging 


1 pc+ioo 


r2(f)r2(f) 

T{ps)r{qs) 


,stl 


( 2 . 22 ) 


If we take the integral contour C_ in fig. 1, then the leading large p contribution comes 
from the residue at s = 0: 


1 


Resr(s)r(-s) 

s=0 


T{ps)T{qs) 


TT—^ -c-^ -2- 

oTT^pq bpq TT^pq 


(2.23) 


In the region Res < 0, the integrand of (2.22) has the following three types of poles: 

(n G Z>o). (2.24) 


2n 2n 

s =-, s =-, 

p q 

The residue at s = — 2n/p is given by 


1 


r2(f)r2(f) 


Res r(s)r(-s) ^^) 2 V 2/ g.M = 


s = —n, 


1 /2n 


dvr^ s=-2nlp 


T{ps)T{qs) 


27rri V n 


CSC 


2TTn\ _2ni± 


P 


r(-2f) 


-e p 


(2.25) 
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The residue at s = —2n/q is obtained by exchanging p and q. The residue at s = —n is 
also computed as 


_ ^ Res ^ ^ J ^ e- 

47r2 s=-n T{ps)T{qs) dvr^n r(—np) T{—nq) 


n/i 


(2.26) 


We conclude that the large p expansion of the classical grand potential takes the form 




4 _ p2 _ ^2 p3 _|_ ^3 


/^+ ^2pg + (2-27) 


r(0) 


do) 


3'K‘^pq 


6pq 


where 


^M2,i(/^) - 


1 /2n 


^M2,Il(/^) - 


n=l 

CO 


27rn V n 


1 /2n 


7^M2, 111 ( 7 '') — ^ 


n=l 

CO 


27rn V n 


CSC 


CSC 


2vrn\ r2(-M) 

e p , 


p ) r(-^) 

2'Kn\ _2nii 


q 7 r(-^) 

(_i)--ir2(-f)r2(-f) _ 


-e 9 


(2.28) 


47r2n r(—np) r(—nt?) 

These expressions reproduce the results in [38]. 

2.4 Semi-classical analysis 

Let us proceed to the quantum corrections to the grand potential. To compute the quantum 
correction to Cp('S)! we use the Wigner transform, as in [2]. For a given operator A, the 
Wigner transform is defined by 


, , /■°° dQ iPoL / Q' 

Av,(Q,P) = J_^ — e « (<?-- 

The trace of A is then given by the phase space integral 


A 


»4) 


(2.29) 


/' 


dQdP 




(2.30) 


(2.31) 


Let us apply the Wigner transform to the inverse of the density operator 

a = = (2™h|)’''(2cosh|)’’(2c„.h|)’''. 

As shown in appendix A, the Wigner transform of O is given by 

0w(Q, P) = ^4cosh^ Y ~ 4sin^ T^kdp '^ ^2cosh^^ . (2.32) 

Using this result, one can easily compute the WKB expansion of 0-w{Q,P) up to any 
order. We would like to compute the spectral trace 


Cp(s)=Trp^ = Tra- 


(2.33) 
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Expanding O ^ around Ow as in [2, 42], the Wigner transform of O ® is computed by 

(2.34) 


(a-^)w = ^ - 0w)nw, 


r=0 


where {s)r is the Pochhammer symbol. To compute the summand in (2.34), one needs 
the Wigner transform of operator products. The Wigner transform of a product of two 
operators is computed by 

{A ■ B)y/ = j4w * By/, (2.35) 

where the Moyal product * is defined by 


A-k B -.= A{x, p) exp 


EE(-i)’ 

71=0 m=0 


ih 


(dxdp — dpdx) 


mhl[2 I ^ P 


B{x,p) 

A{x,p)d^d^-^B{x,p). 


(2.36) 


In this way, one can compute TrO"^ up to any desired order, in principle. However, the 
integral appearing in (2.34) for general s is complicated and hard to evaluate. Practically, 
it is sufficient to compute it for s € Z>o. Here, we use an interesting idea in [38]. The 
quantum correction (pL) can be constructed by acting a non-trivial differential operator 
on the classical one: 


j(")(/r)=p(")jW(/r), (n = l,2,...), (2.37) 

where pl") is a differential operator of p. Its explicit form up to n = 2 was computed 
in [38]. An efficient way to fix this differential operator is as follows. We first compute 
the expansion of J7 ^"'^(k) around k = 0. This can be done by using the formula (2.34) for 
s G Z>o. Taking an ansatz of the form of we try to fix unknown parameters to match 
the first several coefficients of If the ansatz is correct, the obtained result must 

reproduce higher order coefficients. In this way, one can verify the obtained operator up 
to any desired order. Using this method, we have indeed fixed the differential operator up 
to n = 17. The result up to n = 4 is given in appendix B. 

Using this method, we finally find that the semi-classical grand potential has the 
following large p expansion 


77wkb(M) k) = 77pert(Ai, k) + Jm2,i(i^, k) + JM2,ii(li, k) + 77M2,iii(li, k), 
where 77pert(/^;^) is the perturbative grand potential given by 


(2.38) 


with 


T^ert (/^) k^ — 


Cp,,{k) = 


Cp,q{k) 3 


TT^pqk 


p + Bp^q(k)p-\-Ap^q(k) 


Bp,q{k) — 


Qpqk 


\{k), 

(2.39) 

pqk 
^ 24 ■ 

(2.40) 
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The constant part Ap^q{k) is a complicated function of k, whose exact form was conjectured 
in [37] 

= \ {P^^ciqk) + q‘^Ac{pk)) , (2.41) 

where Ac{k) is given by [28] (see also [44]) 

2C(3) 


Ac{k) = 


vr^fc 


1 - 


l3 \ roo 

dx—rz -r log(l — e“^*). 


16 / vr 


okx _ 


(2.42) 


2npi 

q . 


There are three types of exponentially suppressed corrections with the following forms 

CO CO 

_ 271 jJ, » ^ 

an{p,q,k)e p, JM2,ii(/^,fe) = / ^an{q,P,k)e 

n=l 
oo 

JM 2 m{k^, k) = '^ ^nip, q, k)e 


^-nfi 


n=l 


(2.43) 


n=l 


Note that an{p,q,k) is not symmetric in p and q, while (3nip,q,k) is symmetric (namely, 
an{p,q,k) 7 ^ an{q,p,k), I3n{p,q,k) = /3n{q,P,k)). Our task is to fix these coefficients. 

In the later analysis, it is convenient to introduce a function D{s,p,q,k) by 


T>(e*'^) = D{s,p,q,k)e">^, 
where P is a generating function of 

CO 

V = 1 + 

n=l 


(2.44) 


(2.45) 


The definition (2.44) means that D{s,p, q, k) is obtained by replacing 5^ in T) by s. Then, 
the WKB expansion of the spectral zeta function is simply given by 


Cwkb(s) = — D{s,p,q,k), 


(2.46) 


where the classical part Cp^\s) is given by (2.21). Also, the membrane instanton coefficients 
in (2.43) are generically given by 


an{p,q,k) = 


where an\p,q) and /3n^{p,q) are the classical parts in (2.28): 

= SfS ( „ j (—) ivS) ■ <>’• = 4n^n n-np) r(-n9) ’ 

^ (2.48) 

Thus our goal is to determine D{—2n/p,p,q,k) and D{—n,p,q,k). It is useful to notice 
the following symmetric property: 

D{s,p,q,k) = D{-s,-p,-q,k). (2.49) 

One can show this by using p\p^_p q^_q = p~^ and Cp-i(-s) = Cp('S). 


an ^ (P, q) T^{ 2n 


k 


P 


D[ - ,p,q,k], /3nip,q,k) = 


Pn\p,q) 


k 


D{-n,p,q,k), (2.47) 


,(o) 
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2.5 TBA approach 


There is another approach to compute the grand potential by using the so-called TBA 
equations. In [45], the semi-classical expansion of the ABJM Fermi-gas was computed in 
this approach. In the present situation, for p = 1,2, we can use this method. For p = 1, 
the density matrix is given by 


^ ^-lU{xi)-^Uix2) 

P{xi,X 2 )L^ = ^ ^ ’^ith U{x) = qlog 


2 cosh 


X 

2J 


(2.50) 


27rk 2 cosh (' 

The kernel (2.50) has the same form as the one in [46], and one can immediately use the 
result there. The functional equations are given by 


7rik\ 


/?+ ( X -\—— j i?_|_ ( X - — ) exp 


Trik 


Trik\ 


Ulx + —]+Ulx- — ] 


2 J 

Trik\ 


Tiik\ 


2 J 


2 J 

TTik\ 


= l + rf{x), 


Tj \x —— + rj\ X -— = —kR+{x 


2 7 


and 


R-{x + ^) R-{x-^) 

R+{x + ^) ~ R+{x - ^) 


= 2ikdx arctan 77 ( 2 ;). 


The grand potential is computed by 


dnJiK, 


\p=l 


Airk 


/ oo 

dx{R+{x) + R-{x)), 

-OO 


(2.51) 


(2.52) 


(2.53) 


(2.54) 


where an integration constant is fixed by the condition = ^,k) = 0 . 

For p = 2, the density matrix is 

^ XI— 3:2 

p{xi,X2)\p=2 4 ^ 2 ^ ( 2 cosh a)( 2 coshf )®/^sinh2^' 

In this case, the density matrix is different from the one in [46]. Nevertheless, as explained 
in appendix C, we find the following functional equations determining the grand potential 
for p = 2 in a similar way to the p = 1 case, 

^{x + Trik)^{x — Trik) = rf{x) — 1 , 

r]{x + Trik) + r]{x — Trik) = 2^{x) cosh2r(x), 

w(x + iTik) w{x — Trik) 

——^-— - -— = 2idx arccothry(x), 

^[X + TTlk) ^[X — TTlk) 

where ^(x), p{x) and w{x) are unknown functions, and r{x) and t{x) are given by 


(2.55) 


r(x) = arcsinh yM, 

The grand potential is then given by 

i9ki7(k, A;)| 9 = — / dxw{x)sm]i^r(x). 

” ^ '7TK. I 


(2.56) 


(2.57) 


The functional equations (2.51), (2.52) and (2.55) can be solved systematically around 
k = 0. Therefore one can compute the WKB expansion of the grand potential. 
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2.6 Non-perturbative corrections: worldsheet instantons 

So far, we have considered the semi-classical analysis, which is perturbative in the sense of h. 
As explained in [2, 47], the grand potential receives quantum mechanical non-perturbative 
corrections in k. These non-perturbative corrections are caused by the worldsheet instan¬ 
tons in the dual string/M-theory and invisible in the semi-classical analysis.^ In the case 
of ABJM Fermi-gas, fortunately these corrections can be predicted with the help of the 
topological string on local x P^. Interestingly, for some special cases, the Fermi-gas 
system is related to a quantum mechanical system associated with the topological string 
[48] on certain CY, as will be seen later. In these cases, it will be possible to predict the 
worldsheet instanton correction, as in the ABJM Fermi-gas. However, in general, we do 
not know such a connection, and do not have a systematic treatment of these corrections 
so far. One approach to compute them is to consider the matrix model computation in the 
’t Hooft limit, as was performed in [9] for the ABJM matrix model. In appendix D, we 
compute the planar free energy of the (!,<?) model and find the worldsheet instanton effect 
in the planar limit. 

Following the argument in [2, 47], one can estimate an order of such a non-perturbative 
correction. Let us consider the classical Fermi surface with energy E: 


H{P, Q) = p log ( 2 cosh 


P 


+ q log ^2 cosh ^ ] = E. 


(2.58) 


This gives an algebraic curve in the phase space. By rescaling P = pP', Q = qQ', this 
expression becomes 

p log ^2 cosh -|-glog ^2cosh = E. (2.59) 


In the large E limit, we can approximate the Fermi surface as 


H{P', Q') = FI ~ log ( 2 cosh -|- log (2 cosh ^ 


(2.60) 


up to exponentially suppressed correction. This approximated Hamiltonian leads to the 
equations of motion 


■, dH 1 ^ P' ■, dH I , Q' 


dp, 2 2 

On the equi-energy orbit H{P', Q') = E, Q' becomes 


(2.61) 


Q'= iJl - 16cosh^ ^e-2^. (2.62) 

Zi \ z 

^Very recently, a new scenario was proposed in [42]. This scenario states that the non-perturbative 
correction to the grand potential is produced by the perturbative resummation of the spectral zeta function 
via the integral transform (2.18). It would be interesting to explore the worldsheet instanton correction in 
this approach. 
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The solution of this equation of motion is given by 

tanhf = ™sn(i.™). (2.63) 

where sn(u, m) is the Jacobi’s elliptic sine function and = 1 — 16e“^^. As the function 
of t, this has the real period wi and the imaginary period u; 2 , 


wi = 16K(m), UJ2 = 8iK(m'), 


(2.64) 


where is the complete elliptic integral of the first kind and = 1 — m?. Now we 

consider the complexified Fermi surface, in which we regard Q' and P' as complex variables. 
Then, the complexified Fermi surface (2.58) determines a Riemann surface, and we have 
two kinds of periods associated with this Riemann surface [2]. One of them computes the 
volume surrounded by the surface (2.58). We refer to this cycle as the “B-cycle” and to 
the other as the “A-cycle”, following [47].® The large E behaviors of the periods can be 
easily estimated. Noting 

P'{t) = — log ^cosh^ oJi = 2E + (2.65) 

the period along the B-cycle is given by 

p "I 2E 2E O L '2 

(^PdQ = — P'{t)Q'{t)dt + 0{e~l^,e~^) = - + 0{E). (2.66) 

Jb pq J-^i/2 pq 

In order to compute the A-period, it is convenient to use the variables 

Q' = i6, t = IT. (2.67) 


Then we obtain 



i 

pq 



P'{t)0{T)dT + 0{e p 



SwiE 

pq 


+ 0 ( 1 ). 


( 2 . 68 ) 


As explained in [2], the quantum mechanical instanton effect is related to the A-period. 
The leading order of this correction is 



if , „ 


4.E' 

exp 


= exp 

pqk 


(2.69) 


4m 

This means that the grand potential receives the non-perturbative correction of order e . 
We conclude that the worldsheet instanton correction in the present case is expected to 
take the form 

CXD 

_ 4771 . . 

dm{p,q,k)e (2.70) 

m=l 

We observe that this expectation is precisely consistent with the exact computation of the 
partition function for various {p, q, k) and the planar solution of the (1, q) model analyzed in 
appendix D. For the very first few coefficients, we can guess the exact forms of dm{p, q, k), 
as given in the next two sections. In addition to the worldsheet instantons, there also exist 
bound states of the membrane instantons and the worldsheet instantons. Computation of 
these bound state contributions is beyond the scope of this work. 

®Note that this convention is opposite to the one in [2]. 
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2.7 Exact computation of the partition function 

In this subsection we present our algorithm for the exact computation of the partition 
function with fixed integer p, which is a simple generalization of the ABJ(M) case [11, 13, 
17]. First let us recall the formula for the grand partition function 


S(k, k) = exp 


CXD / 


n=l 


n 


(2.71) 


where we define multiplication and trace of two matrices pi , p 2 as 

poo poo 

PiP 2 {xi,X 2 )= dy p{xi,y)p 2 iy,X 2 ), Trpi = / dypi{y,y). 


(2.72) 


This formula tells us that we can exactly compute the canonical partition function with 
the rank N if we find exact values of Tr p"' with n = 1, • • • ,N. Therefore, below we explain 
how to compute the values of Tr p^ exactly for integer p. 


2.7.1 The case of odd p 

When p is odd, the density matrix is given by 


p{xi,X2) 


1 _ 1 _ 1 _ f{xi,X2) 

2{p-l)\TTk (2 cosh (2 cosh 2 cosh 


where 


p-i ^ 
2 


f{xi,X2) = 

j=i L 


3^1 - 372 y (2j - 1)^ 
2Trk 


p-i 

j,j'=o 


Then, we rewrite the density matrix as 


p-i 


p{xi,X2) = 


E{xi)E{x2) 

M{xi) + M{x 2 ) 

.7J =0 


fjj'Xixi, 


where 


E{x) = 


e 2 fc 


(2 cosh 

This relation is also schematically represented by 


M[x) = 2{p — l)\T:kek . 


p-i 


{M,p}= ^ f^y{x=E)®{x^'E). 

jJ'=o 


(2.73) 

(2.74) 

(2.75) 

(2.76) 

(2.77) 


Here we regard p, M and x^E as the symmetric matrix, diagonal matrix and vector, re¬ 
spectively, whose indices are the coordinates {xi,X 2 )- Applying this relation iteratively, we 
find 

n—1 p—1 

{M,p-] = ^^(-1)' /,y(pV^) ® (p"-'-V'e), (2.78) 

^=0 j,j'=o 
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which is equivalent to 


n—1 


P-1 




Here cj)^^\x) satisfies the recursion relation 

^•^ ^ r°° 

</>? (a;) = 


Eix) 7_oo 27rA: 


Pyp{x,y)E{y)(t)^^\{y), 


with the initial condition 


,(3) _ 


= X-'. 


Once we know the series of functions (f)^p{x), we can compute Trp”' systematically. 


2.7.2 The case of even p 

For even p, the density matrix is given by 


p{xi.,X2) 


1 _ 1 _ 1 _ f{xi,X2) 

4(p-l)!7 r2A:2 (2 cosh (2 cosh f 2 sinh 


where 


2-1 r 
2 ' 


f{xi,X2) = {Xl - X2) 

j=i L 

Then, we rewrite the density matrix as 


Xl - X2 

2'Kk 


+r 


p-i 


p-i 

fjfx\x(. 

j,j'=o 


p{xi,X2) = 


E{xi)E{x2) 

M{xi) - M{x 2 ) 

.7J =0 


fjj'Xixi, 


where 


E{x) = 


e2k 


M{x) = 4(p — ly.Ti'^k'^ek . 


(2 cosh 

This relation has the similar structure as in the odd p case: 


p-i 


W,p\= Y fjj'{x^E)®{x^'E). 
iJ'=o 


Hence, a similar argument leads us to 

n—1 


. n—l p—1 

p'^{xuX 2 ) = M{xi)-M{x 2 ) fjj'PP {Xl)pp]i_^{x 2 ) 


where (j)P{x) satisfies formally the same relations (2.80) and (2.81) as in the odd p 


(X2). 

(2.79) 

(2.80) 

(2.81) 

(2.82) 

(2.83) 

(2.84) 

(2.85) 

( 2 . 86 ) 

(2.87) 

case. 
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Figure 2. The free energy {—\ogZ) is plotted against for some {k,p,q). The straight lines 

show the perturbative free energy (— log-Z’pert)- 


2.7.3 Free energy for various {p,q,k) 

Using the method explained in the previous subsections, we have computed the exact values 
of partition functions up to certain N = A^maxj for various {p,q,k). For example, for the 
case of {p, q, k) = (1, 2, 2), (2, 3,1), (2, 3, 2), we have computed the exact partition functions 
up to A^max = 66, 32, 33, respectively. Some examples of the exact values of partition 
functions can be found in appendix F. These exact data are very useful to extract the 
instanton corrections as in [11]. 

From the general argument in the Fermi gas approach [2], in the large N limit the 
partition function Z{N,k) behaves as 

Z{N, k) = Zpert(iV, k)-[l + Oie -'^)], (2.88) 


where the perturbative part is given by the Airy function [2, 10] 

Zpert(A^, fc) = Cp,gik)-^/^e^^’^^’^^Ai\Cp,g{k)-^/^{N - Bp^g{k)) 


(2.89) 


The constants Ap^g{k), Bp^g{k) and Cp^g{k) appearing in Zpe^t{N,k) are none other than 
the coefficients of the perturbative part of the grand potential (2.39). In fig. 2, we plot the 
exact values of the free energy (— log Z) for some {k,p, q) with the perturbative free energy 
(— log Zpej-t). We can easily see that the exact free energy shows a good agreement with the 
perturbative free energy since their difference is exponentially suppressed in the large N 
regime. We also observe that the free energy scales like for large N as expected from 
the AdS/CFT correspondence as found earlier in [2, 49, 50]. The perturbative free energy 
also contains the log-correction | log N in subsubleading large N correction as expected 
from the one-loop analysis on the gravity side [51]. 
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2.8 A comment on the gravity dual 

The (p, (7)-model is expected to be the effective theory of N M2-branes on the orbifold 
{Cy{Zp X Z,))/Zfc [32]: 

27ri _ 27ri 

(t>A : {zi.,Z2.,z^,Zi) ~ (e zi,e Z2,z^,za), 

27r^ 27r'i 

4>b ■ {zi, Z 2 , Z'i, za) ~ {zi,Z 2 ,e p Z 3 ,e p Z 4 ), 

27ri _ 27rz 27r’i _ 27Ti 

4>c ■ {zi,Z 2 ,zs,za) ^ {e'^i Zi,e Z 2 ,e i^p Z 3 ,e >^p za). (2.90) 


This implies that the (p, (7)-model is dual to M-theory on AcISa x (5’^/(Zp x Zg))/Zfc with 
the metric 


ds^ 



2 

AdSi 


+ E?ds 


2 

(S7(ZpXZ,))/Zfc> 


(2.91) 


where 


R = {327:^ kpqN)^/^ Ip. 


(2.92) 


Since this background has many nontrivial 3-cycles, we could have discrete holonomies 
of the 3-form potential along the cycles as in the ABJ theory [5]. For Imamura-Kimura 
type theory with equal ranks of gauge groups (without fractional branes), the discrete 
holonomies depend on the ordering of 5-branes in its type IIB brane construction. This 
has been studied in detail in [52] by analyzing monopole operators in general Imamura- 
Kimura type theory. According to the formula in [52], we expect that the gravity dual of 
the (p, g)-model does not have the discrete holonomies. 

There are some predictions on the free energy — log Z from the gravity side. First the 
free energy of the classical SUGRA with the boundary is given by (see e.g. [53]) 


Fsugra = p.93) 

Also by one-loop analysis of the lid SUGRA on AdSA x Xj with the smooth 7d manifold 
Xj, it is known that the one-loop free energy contains the following universal log-correction® 
[51] ^ 

-log A. (2.94) 

On the CFT side, this behavior comes from the the Airy functional behavior (2.89) in the 
perturbative free energy. 

Next we give some comments on nonperturbative effects. Let us first recall the ABJM 
case. For the ABJM case (p = (; = 1), if we identify the M-theory circle with the orbifolding 
direction by and shrink the circle, then the lid supergravity on AdSA x becomes 

the type IIA supergravity on AdSA x CP^. In the type II superstring on AdSA x CP^, 
we have worldsheet instanton effect, which comes from fundamental string wrapping the 
nontrivial 2-cycle CP^ in CP^ [6]. From the M-theory viewpoint, this corresponds to an 
M2-brane wrapping the non-trivial 3-cycle S^jZ^ in S^/Z^. For the general {p,q) case, we 

® When Xr has fixed points as in our case, there might be extra massless degrees of freedom and the 
logarithmic behavior could change. However, the agreement to the CFT side implies absence of such extra 
contributions. 
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also expect that there are similar non-perturbative effects as in the ABJM case. Note that 
the orbifold (2.90) includes the nontrivial 3-cycle {S^/CZp x Zq)/Zfc), which is obtained 
by taking Z 2 = Z 4 ^ = 0, for example. This implies the presence of non-perturbative effect 
coming from M2-brane wrapping this cycle, whose weight is given by^® 


exp 


-TmsVoI {{S^/{Zp X Zg))/Zk) 


= exp 


2kpq 


exp 



(2.95) 


Since the 3-cycle becomes two-dimensional in the large-A: limit, we expect that this effect 
corresponds to the worldsheet instanton effect described by the fundamental string wrap¬ 
ping a 2-cycle in the type IIA superstring theory. One can see that the weight (2.95) of 
the worldsheet instanton effect computed from the gravity side correctly reproduces the 
weight (2.69) obtained by the matrix model, after changing the variable from the canonical 
to the grand canonical ensemble. This is also consistent with the planar solution of the 
(1, q) model computed in appendix D. 


3 Results on the (1, g)-model 

In this section, we summarize some explicit results in the case of p = 1. In this case, the 
system can be thought of as a one-parameter deformation of the ABJM Fermi-gas by q, 
or the deformation of the Nf matrix model [26, 27] by k. Therefore, the results for k = 1 
(and q = Nf) must reproduce those in the Nf matrix model. Similarly, in the limit g —)• 1 
(with general k), the results must also reproduce those in the ABJM Fermi-gas. 


3.1 Membrane instanton corrections 

Let us first consider the membrane instanton corrections. First of all, one notices that the 
classical membrane instanton corrections j and jjj in (2.28) are divergent in the 
limit p —)■ 1. As was shown in [38], these divergences are, however, canceled by each other. 
One finds that after the cancellation, the finite part is given by 


^0) 
'-'M2,1 




= lim(j7-4'5 1 


+ '^M2,iii) = i9)k + 4°^ {<l))e 

ri=l 


(3.1) 


where 

j(0)(q) = _j_pn\ TH-nq) 

27r2n \n ) T{-2nq) ’ 

= - 4 ;^ Q r{-Z^q) ^ + 2nq{il,{-nq) - ilj{-2nq))]. 

Note that the tension Tm 2 of the M2-brane is given by Tm 2 ~ l/(47r^/p). 


- 18 ~ 















Here Hn is the n-th harmonic number, and ipiz) = 5zlogr(z) is the digamma function. 
There is no limit problem for J^n'- 


^M2,Ii(A^) = 1™^M2,II = ^Oin\<iy 

^ n=l 


2nii 

q 


a® (9) = 


1 2n 


CSC 


2Tm\ r (—^) 


q J Y{-f 


27rn \ n 

Therefore, in the case of p = 1, the large fi expansion (2.38) reduces to 




where 


— ’jj' ^ + Bi^q{k)fl + Ai^q{k) + J7 m 2 ,i(M) k) + i 7 m 2 ,ii(a'') k), 


^M2,l(/i, k) = ^(7n(g, k)n + 6n{q, k))e 2”^, 

n=l 
00 

k) = '^ Oniq, k)e 


2nfj, 

q . 


n=l 


(3.3) 


(3.4) 


(3.5) 


Note that an{q,k) = an{q,ltk), not an{l,q,k). Acting the differential operator in 
appendix B on the classical grand potential, one can find the WKB expansion of each 
coefficient. 


Coefficients of e 

compute 


To find the WKB expansions of 7n((?, k) and 6n{q, k), we need to 

p(n)(^g-2nM) p(n)^g-2n/.)_ (gg) 


The computation of 7n(g, k) is relatively easy, compared to the other instanton coefficients 
an{q, k), 6n{q, k). It takes the form 


'ln{q,k) = 


7n\q) 

k 


D(-2n, l,q,k), 


(3.7) 


where, as mentioned in the previous section, D(—2n, 1, q, k) is obtained by replacing in 
the differential operator T> by —2n. Therefore its WKB expansion is 




+ 


24(-l + 2ng) 

n?{l + 2n)q^ (4 — 24n + 4(7 — 3nq + lAn^q) . 


(3.8) 


5760(-3 + 2ng)(-l + 2ng) 


inky + Oik^). 


Since we have fixed up to n = 17, we can compute the WKB expansion up to 0{k^^). 
All of the following results reproduce the correct WKB expansions up to this order. From 
the WKB data, we find analytic expressions for n = 1,2: 


D{-2, 1, q, k) = 2 F 1 ^ - g; sin^ , 

T>(-4,1, q,k) = ^ 2 F 1 (^-q, -q] ^ - 2q; sin^ ^ 2 F 1 (^-q, -2q; ^ - 2q; sin^ . 

(3.9) 
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One can check that, in the limit /c —^ 1, 7i((?, A:) and 72(9, fc) given by (3.7) with (3.9) 
reduce to 


71(9, 1 ) = 


1 r^(-g/2) 

27r2 T{-q) 


72(9,1) 


_ J_ 7 + £!(zh, 

dvr^ Y cosirqj T{—2q)' 


(3.10) 


which are in perfect agreement with the known results of A^j-matrix model in [28]. 

The constant part dn{q,k) is more involved. Both of (3.6) contribute to 6n{q,k). 
The latter contribution is just D{—2n,l,q,k). A simple computation shows that the p- 
independent contribution of the former is given by dsD{s, 1, q, k)\s=- 2 n- We conclude that 


Sn{q,k) 


{q) 


(^D{-2n, 1, q, k) + dsD{s, 1, q, k)\s=- 2 n^, 


(3.11) 


where the WKB expansion of dsD{s, 1, q, k)\s=- 2 n is given by 


dsD{s, l,q,. 


\s=—2n 


nq^ (—1 — 3n + ng + 4n^g) „ 


, [-3(-3 + 16n + 60n=) + (9 - n + 126n^ ( 3 . 12 ) 

+ 510n^)5 - 2n(8 + 13n + A?,r? + 2l2n^)q^ + 2r?{2 + 5n 

+ 12n^ + 56n^)g^] (vrA:)^ + O(k^). 


It is difficult to guess an exact form of this expansion even for n = 1. 

Coefficients of Next, let us consider the second type correction 77 m 2 ,ii- The 

coefficient an{q,k) is given by 


aniq,k) = 


Oin'’ (g) 
k 


D ( -‘^,l,q,k ). 


(3.13) 


The WKB expansion of D{—2n/q, 1, g, k) is given by 


D ( - —, l,g, A: ) = 1 - 


n^(2n + g) 
24g(-l + 2n) 


(vrA;)^ 


+ 


n 


''{2n + q) (—24n + 14n^ + 4g — 3ng + 4g^) 
5760g2(-3 + 2n)(-l + 2n) 


(3.14) 


{T:kf + 0{k^). 


It is not easy to find out an exact expression of this expansion, but for n = 1 we find a 
surprisingly simple expression in terms of the g-gamma function. 


D 


-1 I a }^ = + Vg) rq(l + 2/g) 

g’ r(l + 2/g) r2(l + l/g)^ 


q = 


(3.15) 


where rq(2;) is the g-gamma function defined in appendix E. As we will see in the next 
section, this conjecture comes from the analysis of this instanton coefficient for the 1/g G Z 
case (4.21). When kq G R, there is a subtlety of the definition of the g-gamma function due 
to jqj = 1. As explained in appendix E, one can define the g-gamma function with jqj = 1 
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by regularizing the infinite product using the zeta-function regulariation. Using the result 
(E.12) in appendix E, one obtains the all-order WKB expansion 


D ( l,q,k] = exp 


^ ( _ ■\\n—l]D 

£ 2n(2n + 1)1 (l + q~^) - 52n+i (l + 2q-^)) {Trkqf^ 

.n=l \ /■ 


Trk 

= -cot 


TTk\ 


exp 


(-l)-iB 


^ 2n(2n-h 1)! 

.n=l 


{2B2n+i {q - B2n+1 {2q ^)) {nkq) 


2 n 


(3.16) 

where B^ and Bn(x) are the Bernoulli number and the Bernoulli polynomial, respectively. 
Erom the first line to the second line in (3.16), we have used the identity 


Bm(l -h rc) = Bm{w) mw 


rn —1 


(3.17) 


By using an integral representation of the Bernoulli polynomial [54] 

sin 2t:w 


f-C 

B2n+i{w) = {-lT-\2n + l) / 

Jo 


dt ■ 


-t 


2 n 


Iq cosh 27rt — cos 27rrc 
one can perform the resummation of (3.16). For q >2, one easily hnds 

D = ^cot exp[/i(g,fc)], g > 2, 


0 < w < 1, 


where 


h{q,k) = f 
Jo 


dt- 


4sin2^sin^(cosh27rt-h2cos2 2) ^ /sinh(7rA:gt/2) \ 


■log 


0 (cosh 27rt — cos ^) (cosh 27rf — cos ^) ” V T^kqt/2 J 


(3.18) 


(3.19) 


(3.20) 


For 1 < g < 2, to use the integral representation (3.18), we further need to shift the 
argument by using (3.17) with w = 2q~^ — 1. Then we find 


D 




sm^{2-q) 

fi‘2-q) 


exjp[Ii{q,k)], 


l<q<2, (3.21) 


where Ii{q, k) is the same as above. In a similar way, one also obtains the integral repre¬ 
sentation for 0 < g < 1. These integral representations are very useful to understand the 
pole structure of ai{q, k). Since the integrand of /i(g, k) is exponentially dumped for large 
t and does not have singularities in the integral domain for g 7^ 1,2 and hnite k, Ii{q,k) 
takes in a finite real value and hence does not have any singularities for A: > 0. Thus the 
singularities of C(i{q, k) come only from the cotangent factor in (3.19) or (3.21). In figure 3, 
we illustrate these for q = 3. 


3.2 Worldsheet instanton corrections 

In the case of p = 1, the worldsheet instanton corrections take the form 

00 

_ 4772 

dmil,q,k)e 9'='^. (3.22) 

771=1 
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/i(3,k) 



simple poles at even k. 


From a consistency with the results for various integral (g, k) (see appendix F), we conjec¬ 
ture the exact form of worldsheet instanton coefficients for m = 1,2: 


di{l,q,k) 


^2(1,9, fc) 




(3.23) 


Note that these are also consistent with the planar free energy computed in appendix D. 
In the limit q ^ 1, these precisely reproduce the worldsheet instanton corrections in the 
ABJM Fermi-gas computed in [11]. Also, for q = 2, (3.23) reproduces the worldsheet 
instanton coefficients of {p,q) = (1,2) model found in [38]. 


3.3 Pole cancellations 

Since we have determined some of the instanton coefficients analytically, we can see the 
pole cancellations in some special limits of {q, k) beyond the semi-classical approximation. 
These are important non-trivial tests of our conjectures. 


ABJM limit. In the limit 1, all of an(q, k), 7n(^, k) and dn{q, k) are divergent. Let 
us see the cancellation of these divergences for n = 1. We first notice that Ii{q, k) behaves 


as 


\iT[vh{q,k) =0{{q-lf). (3.24) 

q^l 

Using the integral representation (3.21), one finds that the divergence of ai{q,k) is given 
by 


lim ai{q,k)e 

q^l+ 


2/i 

q 



2{q-iy 


+ 


vr/c TTk\ 1 

u + I -cot — - 

^ 4 2 J q-1 


+ 


Trk Tik 

1-cot — 

2 2 


1 7r^(4 — k‘^) 

^+2+ 16 


Trk Trk 

— cot — + 0(g- 1) 




(3.25) 


The divergence of 71 (g, k) is 

lim 7i(g,/c) 
q^l 


4 




cos 



-^ + l + 0{q-l) 


(3.26) 
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Thus we get 


lim (ai{q,k)e “ + ii'yi{q,k)e 
g^l+ V / 


4 / 7rk\ \ 1 

VT ) Uiq-iy 


, nk Trk 


q-1 


=-2m 


+ {ai{k)fi^ + bi{k)fi + ci{k) + 5ci{k))e + 0{q - 1), 


(3.27) 


where 


ai{k) = 
ci{k) = 


4 irk 


2 Trk Trk 

bi{k) = - cos —cot—, 

IT 2i 2i 


2 5A: 1 Tik k 2 '^k\ irk 

- H- 1 — cot- 1 — CSC — cos —, 

3A: 12 vr 2 2 2 / 2 ’ 


(3.28) 


and 


6ci{k) 



(I , 7r^(4-fc^) 
V2 16 



ci(fe). 


(3.29) 


The divergence of the /r-dependent part is precisely canceled. Furthermore, the coefficients 
ai{k) and bi{k) in the finite part perfectly agree with the resnlts in the ABJM Fermi-gas 
[11]. The divergence of the //-independent part must be canceled by 5i{q,k). This means 
that 6i{q,k) must behave as 


lim 5i{q,k) 

q^l+ 



1 f Trk Trk\ 1 

L2(g-i)2 + (^"T"°*Tj 7 ^ 


6ci{k). (3.30) 


This is regarded as the constraint for 6i{q,k). 


/c —)• 2 limit. Let us consider another limit A: —)• 2. In this limit, the coefficient ai{q,k) 

_ 2/j, 

of e ”7 diverges, as shown in figure 3. This divergence must be canceled by the leading 
worldsheet instanton correction of order e . It is easy to find 




lim di(l, q, k)e 9 * = 
k^2 


2q 


TT sin - k — 2 
<? 
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2/t -I- g + TT cot ^ 


vr sm ■ 


+ Oik-2) 


27r 

e 9 . 


(3.31) 


Using the integral expression (3.19) or (3.21), we also find 


_ 2 //. 

\im ai{q,k)e ~ = 
k —^2 


2q 


cot ■ 


TT sin - k — 2 
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sm - 
<? 


+ Oik-2) 


e 9 


(3.32) 


The singular parts are indeed canceled as expected. The finite part is finally given by 


lim 

k^2 


( _2ii _4 m\ 

[ai{q,k)e 9 +di{l,q,k)e 9fej 


2u, + q _ 22 i 

- e 9 . 

TTsin - 


(3.33) 


This correctly reprodnces the coefficient (F.13) of e for the {p,q,k) = (l,g, 2) case 
(see appendix F). 
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fe —>• 2n limit. More generally, ai{q,k) has a pole at k = 2n (even integer). From the 
integral expression (3.19), we find that ai{q, k) behaves in the limit A; —>• 2n as 


lim ai{q, k) = 
k^2n 


2q 


A: — 2n TT sin Y 


n 

1=1 


sm 


T<-{n+j) 

nq 


sin ^ 


+ 0 ( 1 ). 


nq 


This pole should be canceled by the worldsheet ra-instanton 

/ 2^ \ 

lim (ai{q,k)e~~i~ + dn{l,q,k)e~~i^ ] = finite. 

fc->-2n V / 


(3.34) 


(3.35) 


One can see that ^2(1; Q, k) in (3.23) indeed has the correct pole structure at A: = 4 satisfying 
the condition (3.35). For n > 3, this pole cancellation condition (3.35) gives the constraint 
for a possible form of dn{l, q, k). 


4 Results on the (p, g)-model 

In this section, we give explicit results for the general (p, g)-model. The basic strategy is 
the same as in the case of p = 1. We compute the WKB expansion of each membrane 
instanton coefficient, and then conjecture its analytic form. Since the WKB expansions 
become much more complicated than those for the (l,g)-model, it is harder to determine 
their analytic forms. To fix the worldsheet instanton corrections, we use the exact results 
for various integral {p, q, k) in appendix F. 


4.1 Exact partition function for N = 2 


We first compute the spectral zeta function at s = 2, exactly. It is easy to find that Cp(l) 
is exactly given by 

1 r2(|)F2(|) 


r 

Cp(i) = / 

J —( 


dx p(x, x) = 


Air^k F(p)r(g) ' 


(4.1) 


Also, Cp(2) can be computed as follows. 


Cp(2) = J dx {a 
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1 


1 
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(2 cosh y )P (2 cosh ^)i (2 cosh ^)p (2 cosh 

f 1 11 

= / dxdy (x| -- \y) {y\ - - \x) 


(4.2) 


(2cosh-|^)P (2cosh-^)P (2cosh 2 )*^ (2cosh ^)'? 


By using the fact that (x| G{P) \ y) depends only on x — y and shifting the integral variable 
x ^ X + y, we find 


Cp(2) = 


dxdy 
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(27rA:)^ 
dx 


{x - y\ 
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(2 cosh 


| 0 ) 
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{x\ 
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|0) P{q,x), 


(2 cosh 1)9 (2 cosh |)9 


(4.3) 


(2 cosh y)^ 
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where 


r 

F{q,x) = / 

J —C 


dy 


_1 1 r(g)^ 

/_oo 2'Kk (2cosh (2 cosh |)9 2Trk r(2Q') 

Using the Fourier transform 


9 9.1 
2 ’ 2 ’ 2 


n; oiTT + 9 ;-sinh - . 


(xl 


(2 cosh ^)P 




p IX P IX 


27rk 2 2'Kk 


we finally find 




|,|;^+9;-sinh^ vrfcx 


9 X 
2 


(4.4) 

(4.5) 


(4.6) 


where we have rescaled the integration variable x —>• 2Kkx. Note that this is an exact 
expression. Recalling the relation (2.15), the exact partition functions for = 1,2 are 
given by 

Z(l, k) = Cp(l), ^(2, k) = -iCp(2) + ^Cp(l)'. (4.7) 

In the case of p = g = 1, Z(2, k) correctly reduces to the exact partition function for N = 2 
in the ABJM theory, computed in [55]. 

For the spectral trace Cp(^) with general i € Z>o, we conjecture that it has a simple 
integral representation 


Ul)= [ ^^^U(g)U(P)e(-^)'*^‘^«^^U(g)U(P)e(-^)' '^^Q^PV{Q)U{P) ■ ■ ■ e^^^‘^^PV{Q)U{P). 
I zTTri 

(4.8) 


Here, for simplicity, we have introduced the notation V(Q) = (2 cosh ^)“'^, U(P) = (2 cosh y )“^, 
and the derivatives in (4.8) act on all functions on their right. One can show that, for i = 2, 

(4.8) indeed agrees with (4.6). Although we do not have a proof of (4.8) for £ > 3, we have 
checked that this conjectured expression (4.8) correctly reproduces the WKB expansion. 


4.2 Membrane instanton corrections 

In this subsection, we consider the membrane instanton corrections. 


Coefficients of e Let us first consider the coefficient of e in (2.43). The WKB 
expansion of D{—n,p,q,k) up to 0{k‘^) is given by 


D{-n,p,q,k) = 1 - 


n^(n^ — l)p'^q^ 


{Kky 


+ 


96( —1 + np)(—1 + nq) 
n^(n^ — l)p^g^ (SOn — 24p — 24n^p — 24g — 24n^g + llnpq + 7n^pq) 
92I60(—3 + np)(—I + np){—3 + nq){—l + nq) 


{Kk)‘^+ 0{k^). 

(4.9) 

One notices that D{—l,p, q, k) does not receive the quantum corrections: Il(—l,p, g, k) = 1. 
This is indeed the case. Since Cp(l) does not receive any corrections, we have 11(1,p, g, k) = 
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1 for any {p,q). Using the reflection symmetry (2.49), we conclude that D{—l,p,q,k) = 1. 
Therefore I3i{p, q, k) is exactly given by 


p{^\p,q) _ 1 r^(-2)r^(-f) 


Pi{p,q,k) = 


k Air^k r{-p) r{-q) ■ 

Moreover, the WKB expansion for n = 2 has the following remarkable structure: 


(4.10) 


Trk\ 


L>(-2,p,g,A:) = ^—^/„(p)/„(g) ( — J , 

n=0 I 1’ k / 

where a generating function of fn{p) is given by 

p pi 


(4.11) 




n=0 


(2n)! 


-;- n: sm z 

2 2 2 ^ 


(4.12) 


More explicitly, fn{p) is given by 




m-j 


(4.13) 


In the previous subsection, we have already computed Cp(2)- From (4.6), one hnds that 
Di2,p,q,k) = 2-^1 ^|, |; ^ + g; - sinh^ (4.14) 

Using the symmetry (2.49), one immediately obtains 


D{-2,p, q, k) = D{2, -p, -q, k) 


r(- 2 p) r 

27rr2(-p) 
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P 


dx B^ ( —-+ix,—- — ix] 2^1 
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1 

V 2 ’ 

2’ 

2 


q-, — sinh^ vr/cx 


(4.15) 

One can understand the factorized structure (4.11) of D(—2,p, q, k) from the expression of 
Cp(2) in (4.8). We emphasize that this expression is valid for any {p,q). In particular, for 
p = 1, 2,3,4, we find the following analytic expressions 
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D{-2, 2,q,k) = - + - 2 F 1 --q; sin^ Trk , 
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D{-2,3,q,k) = ^ 2 Fi( ^ - g; sin^ ^ ^ 2 ^^! 


Trk\ 1 
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q; sin^ 1 ^ 2 -F 1 
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q q 1 2 3vrA: 

g,„ 

q q 1 


, , q: sm 

2 2 2 ^ 2 


2> 2 ’ 2 F,sm 2TTk] . 


(4.16) 

Similarly, if either D{n,p, q, k) or D{—n,p, q, k) (n G Z>o) is known, one can know the other 
by the symmetry (2.49) and the analytic continuation [p, q) (—p, —q)- By matching the 
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WKB data, we find the analytic form of D{—n,p, q, k) for some other cases 


q I 3q . 2 


D{-3, 2, q, k) 2 F 1 ^-q, ^ “ y! sin^ 

D(-3,4, q, k) 2 F 1 y, -9; \ - y; sin^ + y s-F’i ^ “ Y’ 

+ 77 ^^^r^’" 2 ’ 2 "T’®“ 2 vrfcJ+-, 

Z)(-4, 2, A;) 2 F 1 (^-3q, -q] ^ - 2q-, sin^ ^ 2 Y ^-2g, - 5 ; ^ - 2g; sin^ vrA; 


9 


1 


H- 2 Y ( —< 7 , — 9 ;-2(7; sin^ Trk 1 H-. 

35 V 2 y 35 

(4.17) 

Note that when q is an integer, the hyper geometric series in (4.16) and (4.17) are truncated 
to a finite sum, and they are reduced to some combinations of trignometric functions. 

Coefficients of The WKB expansion of the coefficient of up to 0(A:^) 

is 


D 


- — ,p,q,k ) = 1 - 

p 


n^q^{2n — p){2n + p) 
24p(—1 + 2n){—p + 2nq) 
n^q^{2n — p){2n + p) 


[i^ky 


5760p2(-3 + 2n)(-l + 2n)(-3p + 2nq){-p + 2nq) ^ 

- I2p^ - ASn^q + 28n^q - I2p^q + nnp^q){-Kkf + 0{k^) 


(4.18) 


Here we focus on the n = 1 case, i.e., the coefficient of e For the case of q/p = m G 

Z>o, we find that the following expression correctly reproduces the WKB expansion up to 

(4.19) 


One can easily show that this can be rewritten as a contour integral 
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1 + zq 




1 -I mp 


_n=o 4 + zq 




(4.20) 


where q = In particular, when {p,q) = (1/m, 1), this integral can be evaluated 

explicitly by expanding the integrand using the gr-binomial formula. By picking up the 
coefficient of z™, we find 


D 



(2m)! nr.iSinks' 


(4.21) 
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By noticing that (4.21) can be written as a combination of g-factorials, we have arrived at 
the conjecture in the previous section that the coefficient of is given by the g-gamma 

function (3.15) for the general {p, 1) case. Also, for the {p,q) = {2jm,2) case, the integral 
(4.20) can be evaluated exactly thanks to the formula (3.9.1) in [56] 
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m\ 


|2 


(2m)! 


2 <(>i(q 


m! 
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— ,q—,q;q,q-+i; 


nf=i sin2 ^ 


(2m)! n' sin2 ^ 07^' sin^ ^ ^ 


(4.22) 


where 2 '^i(q“) qj z) denotes the basic g-hypergeometric series. This suggests that the 

the coefficient of for the general (p, 2) case is given by a certain analytic continuation 

of (4.22) to a non-integer m. However, compared to the g-gamma function appearing in 
the (p, 1) case, the precise definition of the g-hypergeometric series with |q| = 1 is much 
more subtle (see [57] for some proposal). We leave it as an interesting future problem. For 
g > 2, the integral (4.20) is hard to evaluate explicitly. 

When m = g/p is not integer, the contour integral representation (4.20) is no longer 
correct due to the branch cut of Instead, we conjecture that the coefficient of 

g-2At/p for the general (p, g) case is given by the following integral along the real segment 
0 < z < 1 


D 



^r(-2g/p) dz ( qPipA.) M 

TH-q/p) Jo 


(4.23) 


with 

Fip,k,z) = 2f; (1 + -) Lii_ 2 n(-^)(vrM'"- (4-24) 

^(2n + l)! V2 pj 

Here Li„(z) denotes the polylogarithm. We do not have a proof of this conjecture, but we 
have confirmed that it correctly reproduces the WKB expansion up to 0{k^y. For p > 2, 
one can rewrite F{p, k, z) as the following integral form by using (3.18), 


r-co sin — 

F{p,k,z)= / dt - - - + 2zcos\iT:kpt + z^), p>2. (4.25) 

Jo cosh 2'Kt + cos ^ 

For p < 2, one has to shift the argument of the Bernoulli polynomial in (4.24) to use the 
integral representation. 


4.3 Worldsheet instanton corrections 


As in the same way in the previous section, we can compute the worldsheet instanton 
corrections for various integral (p, g,/c) (see appendix F). From these data, we conjecture 
that the leading worldsheet instanton correction in (2.70) is given by 


di{p,q,k) 


pq 


sin H sin I • 


(4.26) 
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To guess the higher order corrections is not easy. For the (1, g)-model, we conjectured 
the 2-instanton correction in (3.23). We also conjecture the 2-instanton correction for the 


( 2 , q)-model as 


d2{2,q,k) = 


2 


sin 


2 2n 
qk 


.2 TT 


+ 


g(2 + cos g) 


(4.27) 


One can check that in the limit 5 —)• 2 this reproduces the result of {p, q) = ( 2 , 2 ) model in 
[38]. 


4.4 Pole cancellations 

The worldsheet 1-instanton coefficient (4.26) for the general (p, q) case has a pole at k = 2/q 


_ 4/i 

lira di{p,q,k)e pi'^ 

q 


TTSm 




, cot — + q{2p + p) 

k - 7; P 


e p . 


(4.28) 


This pole should be canceled by the membrane 1-instanton coefficient ai{p,q,k) of 
given by (4.23). We have checked numerically that the membrane 1-instanton coefficient 
(4.23) has the correct pole and residue to cancel the pole of worldsheet instanton (4.28). 
Similarly, the pole of the membrane 1-instanton at k = A/q should be canceled by the 
worldsheet 2 -instanton 


lim 

9 


d 2 {p,q,k)e~pi^ +ai{p,q,k)e~ p 


finite. 


(4.29) 


This gives the constraint for a possible form of d 2 {p,q,k). However, it is difficult to nu¬ 
merically calculate the integral (4.23) in the regime k > 2/q, and hence we are unable to 
determine the residue of ai{p,q,k) at k = 4/q so far. Also, it is not clear whether the 
residue of ai{p, q, k) at k = 4/q for the q > 3 case is simply given by trigonometric func¬ 
tions. It would be interesting to find the exact form of worldsheet 2-instanton coefficient 
d 2 {p, q, k) for the general (p, q) case. We leave it as a future problem. 


5 More results in special cases 

In this section, we discuss results for some specific values of {p,q). In some special cases, 
there is a direct connection to the topological strings on certain Calabi-Yau three-fold. 

5.1 Relation to the topological strings 

5.1.1 The (2, 2)-model and the local del Pezzo 

In [39], it was observed that the worldsheet instanton correction in the (2, 2)-model can be 
reproduced by the topological string on the local ZI 5 del Pezzo surface. Here we show that 
the Fermi surface (2.58) with p = q = 2 is indeed equivalent to the mirror curve for the II 5 
del Pezzo. Let us rewrite (2.58) as 

(eW2 + g-P/2)p(gQ/2 ^ ^ 
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For p = q = 2, this reduces to 


+ + Q 26*^ + 4 + 2e ^ + e + e ^ + e ^ ^ . 


(5.2) 


Looking at fig. 1 in [58], one finds that this Fermi surface is identical to the mirror curve 
for the del Pezzo.^^ Following the formulation in [48, 59], the (quantized) mirror curve 
is enough to compute the free energy on the corresponding geometry in the topological 
string. However, one should be careful about the prescription of the quantization of the 
mirror curve. As in [48], a natural way to qunatize the mirror curve is Weyl’s prescription: 

^rP+sQ ^rP+sQ^ ( 5 ^ 3 ^ 


On the other hand, the quantization of the Fermi surface (5.1) leads to 

(e^/2 + g-P/2)p(gQ/2 ^ g-Q/2)q 1^^ ^ gE ^ (5 4) 


where [V') is an eigenstate in the quantum spectral problem. This quantization induces an 
additional factor: 

grPgSQ ^ g-rs«fcgrP+sg^ (5 5) 

Such fc-dependent factors should be taken into account appropriately when one computes 
the membrane instanton correction from the topological string free energy.In particular, 
one should carefully indentify the moduli in the topological string and the parameters in 
the Fermi-gas system. See [47] in more detail in the ABJM case. 


5.1.2 The (1,—l)-model and the resolved conifold 

For {p,q) = (1, —1),^^ the function D{s, 1, —1, k) becomes remarkably simple: 


D{s, 1,-1, k) 


nks 

4 sin ^ 


Then we can compute the WKB grand potential by 


JwKBip, k) 
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ds 

27ri s sin ^ sin^ ^ ’ 


(5.6) 


(5.7) 


where we have used F( 2 :)r(— z) = —^/(z sin (ttz)). By taking the integral contour C+ in 
fig. 1 and picking the poles at s = 4re/A: and s = 2n with n € Z, we find 






4 ' n sin 

n=l 


2 2 7m 
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47r sm 

n=l 


nkn 
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nk nkn 

2nn — 1 -cot- 

2 2 




(5.8) 


^^The local del Pezzo surface corresponds to the polyhedron 15 in fig. 1 in [58]. 

^^Note that the worldsheet instanton correction is determined by the classical mirror curve. We do not 
need the quantization of the curve in the computation. 

^®This case cannot be defined in the original setup. We define this case as a naive analytic continuation 
in the ideal Fermi gas system. 
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The first term is similar to the WS instanton effects while the second terms is similar to 
the membrane instantons. Although each term has poles for rational values of k, these 
poles are actually canceled and the result is finite. 

We can understand this expression from the rehned topological string on the resolved 
conihold as follows. First let us note that the classical Fermi surface for {p,q) = (1,-1) is 
determined by 

p p 


62 + e 2 _ ^ 
~Q 'Q ~ ® 


(5.9) 


6 2+6 2 


1 + 6 *^ + — 6 


P' 


,-2EO'-P' _ 


= 0 , 


(5.10) 


where = {Q, This equation is the same as the mirror curve of the resolved 

conifold (see [60] for instance) and hence we expect that the {p, (;)-model for {p, q) = ( 1 , — 1 ) 
is described by the topological string on the resolved conifold. 

Let us explicitly test our expectation (see also [42]). The free energy of the refined 
topological string on the resolved conifold is given by [60] 


Fieue2;Q) = -f2 




— q _ f n/2'j ’ 

where 

q = e^\ t = 6"^L 

Then the Nekrasov-Shatashvili limit [61] becomes 
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If we identify the parameters as 


61 = vrifc, Q = e + T = ^, -^ = p 


then we find 
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Also, in the standard topological string limit 6 i = —62 with the identifications 

47ri 4u , 

61 = ^, g = e^+-, 


the free energy becomes 


^ , 1 ^(- 1 )’^ 6 ^ 
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By comparing the grand potential (5.8) with (5.15) and (5.17), we easily see 


(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 

(5.17) 


JwKB{k,k) = Ftop{k,p) + 


AFNs(A,r) 


(5.18) 


This structure is the same as the relation between the ABJ(M) theory and refined topo¬ 
logical string on local x P^ [3, 13, 14]. 
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5.1.3 A comment on general case 

For general values of {p, q), we do not find the correspondence to the topological strings on 
the known Calabi-Yau geometries. As in section 2, the WKB expansions of the membrane 
instanton corrections can be computed for any (p, q) even though we do not know the 
topological string counterpart. Let us give a comment how to compute the worldsheet 
instanton corrections systematically for generic {p,q). As seen above, the Fermi surface is 
closely related to the mirror curve of the corresponding topological string. This suggests us 
to regard the Fermi surface (5.1) for general {p,q) as a “mirror curve” of an unconvetional 
geometry. Using the formulation in [62], one can, in principle, compute the genus g free 
energy for this “mirror curve.” It is natural to expect that this free energy just gives the 
worldsheet instanton corrections in the Fermi-gas system. The important point is that the 
formulation in [62] can be applied for any spectral curve even if its geometrical meaning is 
unclear. In practice, however, it is not easy to compute the higher genus correction in this 
way. It would be interesting to test this expectation explicitly. 


5.2 Exact partition function for the (2, 2)-model at A: = 1 


The (2, 2)-model was studied in [39] in detail. Here we point out that the grand potential 
at A: = 1 is exactly related to the topological string free energy on local x P^.^^ The 
modified grand potential for the ( 2 , 2)-model at A: = 1 is given by (see [39] for detail) 

,3 a.,2^a,.^a r 26p‘^ + p + 9/2 


<^2,2 (Z^, 1) — 




2C(3) V + 4^ + 4 

O I r\ ^ I 


TT^ 


TT^ 


+ 2 


TT^ 


e + 


(5.19) 


One notices that this large p expansion is very similar to the one in the ABJM Fermi-gas 
at A: = 2 [11]. In [15], the ABJ(M) grand potential at A; = 1,2 is exactly written in terms 
of the topological string free energy. Recalling this fact, one easily finds that the modified 
grand potential (5.19) is written as 

J 2 , 2 (/i, 1) = ^ [Poit) - tF'it) + + F,{t) + F^^it). (5.20) 


Several definitions are in order. The functions Fb(A) and Fi(t) are the standard genus 
zero and genus one free energies on local P^ x P^, respectively.^^ These are computed in 
a standard way of the special geometry. The function F^^{t) is the first correction to 
the rehned topological string free energy in the Nekrasov-Shatashvili limit. The Kahler 
modulus t is related to the complex modulus z by the mirror map 


t = - log z + 4:Z 4 F 3 


3 3 

l,l,-,-; 2 , 2 , 2 ; 
’ ’22 



(5.21) 


In the present case, the complex modulus 2 is related to the chemical potential or fugacity 
by 


z = 


e->^ = 


1 

K 


(5.22) 


^"^As already seen, the (2, 2)-model is related to the topological string on the local D 5 del Pezzo surface. 
The relation here is probably accidental. 

Since there are two Kahler moduli ti and t 2 in local x P^, the free energy is in general a function of 
these two parameters (ti,t 2 ). Here we denote the free energy in the diagonal slice by Fg{t,t) = Fg{t). 
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(5.23) 


As in [15], the genus zero free energy is written in the closed form 


F''{t) = vr 


K(1 + 16z) 


K(-16z) 

At the large radius point {t —)■ oo), this leads to 




9 


328 


+ TTl. 


777 


Fo{t) = ^ - 2C(3) + 4e-^ - + 0{e-^^) 


(5.24) 


where we have fixed integration constants properly by following [15]. Eliminating t by 
(5.21), one easily finds 

Fo{t) - tF^it) + ^Fo{t) = -^ log^ ^ - 2C(3) + 4(log^ z - log z + l)z 

+ f—261og^ z + logz “ ^ j + ^(8281og^ z + 228 log z + 77 ) 2 ;^ + O(z^). (5.25) 


The free energies Fi(t) and T’^®(t) are also exactly given by 

= -^log[64z(l + 16z)] - ^log 
logz ^2 10 2 I 224 3 


K(-16z) 


TT 

.,4\ 


+ -z -H- z-^ + O(z^), 

12 3 3 9 ^ ^ 


= rf:-logz-^ log(l + 16z) 


1 


12 
log z 


24 


(5.26) 


(5.27) 


2 16 2 512 3 4 , 

-- z-l - 2 ^-+ O(z^). 

12 3 3 9 ^ ^ 

Plugging these results into (5.20), one can check that the large // expansion (5.19) is 
correctly reproduced. 

Exact grand partition function. Once the modified grand potential is known, one can 
compute the grand partition function. The grand partition function is constructed from 
the modihed grand potential by 


uGZ 


J(/i+27rm,fc) 


Plugging the result (5.20) into the summand in this equation, one finds 


^J2,2(/J'+27rin,l) _ gJ2,2(A‘,l) 


exp 


irin T + 27rin £- 

' ^ 3 


(5.28) 


(5.29) 


where 


r = -F^t), £ = ^{tF''{t) - F'{t)). (5.30) 

TT TT 

We have used the identity exp[—= exp[—^1^] with n G Z. Therefore the exact 
grand partition function is expressed in terms of the Jacobi theta function 




(5.31) 
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This expression is useful in // —)• oo. Now we want to analytically continue it to the regime 
/r —>• —oo (or K —)• 0). To do so, we write ^ 2,2 in terms of periods around the orbifold point 
[9, 15]: 


Svr 


1 1 1 


K 


3-^2 > 


2 2 2 


d^ToiX) = 2 ’ 2>2 


2 ' 16 
-„2 ^ 

' + 


K 

16 


n'^ik 


3-C2 


1 1 1 ^ 3 
2’2’2’'^’2’”16 


(5.32) 


where k = k? and is the Meijer G-function. Along the computation in [15], one finds 
that the grand partition function is finally given by 


^ 2 , 2 ( 1 ^, 1 ) = exp 


TT^ 


A%2. 


J-Q - XdxTo + —diTo ]+Ti + Fi 


iNS 




(5.33) 


where 


and 




1 


Fi = -logr/( 2 r) - 2 log 2 . 


(5.34) 


(5.35) 


Now, we can expand the grand partition function around k, = k‘^ = 0. Using the expasions 


TT^ 


- — ( To - A5 aTo + '—dlFo ) = ^ - 


2 


+ 


199k^ 


+ » 


K 


87r2 3847r2 22118407r2 

+ OiK^), 


+ 0(k^ 


32768 524288 

, „ , - 1 , K 1536 - 2567r2 + Dvr^ 2 

logfl2K,f) = 8log« + ^--K 


+ 


368640 + 6144007r2 - 1282887r^ + 20257r6 3 


10616832007r6 


k-^FO{k^), 


we finally get 


(5.36) 


„ , . K 15 — TT^ 2 855 + 757r2 — 167r^ 3 ^ 2 -,/ 4 s 

+ 518400.» 

This precisely reproduces the exact values of the partition function in [39]. 


(5.37) 


6 Conclusions 

In this paper we have studied the partition function of the (p, g)-model on 5^ and in¬ 
vestigated its large N instanton effects by using the Fermi-gas approach. Based on the 
systematic semi-classical WKB analysis, we have found the analytic results on the mem¬ 
brane instanton corrections. The membrane instanton coefficient of the type is related 
to the spectral zeta function Cp(n) by the reflection symmetry (2.49). From the explicit 
forms of Cp(l) ™ (4.1) and Cp( 2 ) in (4.6), we know the exact expressions of the coefficient 
of for n = 1 , 2 . As shown in (4.16) and (4.17), when p is an integer with generic g. 
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the coefficients of reduce to some combinations of the hypergeometric functions. The 
membrane instanton of the type (or is more involved. We found an inte¬ 
gral representation (4.23) of the coefficient of 1-instanton for generic {p,q,k). Very 

surprisingly, for the special case of {p,q) = (l,g), the coefficient of is given by the 

q-gamma functions (3.15). We emphasize that this is quite different from the Gopakumar- 
Vafa type formula [29] in topological string, where only trigonometric functions of h or 1/h 
appear. It is very interesting to understand the physical meaning of this finding better. 
From the observation of the special case (4.22), we speculate that for the general {p,q) 
case, the coefficient of in (4.23) is related to g-hypergeometric functions. 

We have also found some exact results on worldsheet instanton corrections, which 
appear as the quantum mechanical non-perturbative corrections in the Fermi gas, from the 
exact computation of the partition functions at hnite N. We have found the worldsheet 
1-instanton for the general (p, q) case in (4.26), and the worldsheet 2-instanton for the (1, q) 
and (2,g) cases in (3.23) and (4.27), respectively. It would be interesting to understand 
more general structure of the worldsheet instanton corrections for the general (p, q) case. 

We have seen that the apparent poles at the various integral (or rational) values of 
(p, q, k) are actually canceled out between the worldsheet instantons and membrane in- 
stantons, as required. In particular, for the (p, g) = (1,1) case, after the pole cancellation 
the remaining finite part reproduces the known results of the ABJM theory in the highly 
non-trivial way. It is interesting that the quadratic polynomial of p in front of for 
the membrane instanton of ABJM theory correctly appears from the (p, g)-model after the 
pole cancellation. However, this is very mysterious from the viewpoint of bound states. In 
the case of the ABJM theory, one can remove the effects of the bound states by introducing 
the effective chemical potential peff, which is determined by the coefficients of in the 
membrane instantons. However, before the pole cancellation there is no p^ term in the 
membrane instantons. Therefore it seems that there is no natural way to introduce peg in 
the (p, g)-model for generic (p, g). It would be very interesting to study the structure of 
the bound states in the (p, g)-model. 
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A Computing the Wigner transform 

In this appendix, we derive (2.32). The computation is almost the same as the one in [42]. 
By definition, the Wigner transform of O is given by 


r 

Ow{Q,P) = / 

J — c 


dQ' iPQ' ( Q + Q'/2\ 


h 


-e ft 2 cosh 


x(«-f 






) (= 


2 cosh — 
2 


(A.l) 


The last part is written as 


Q- 


Q' 


2 cosh 


P 


Q + 


Q' 


dP' iP'Q 


e ft 2 cosh 


>i\p 


Therefore, we find 


Ow{Q,P)= / sn e 2- 

J — CO 


27r 


( 4 cosh^ ^ + 4 sinh^ ( 2 cosh ^ ] , (A.3) 


(A.2) 


(27r)2 




P' 


/\ p 


where we have rescaled the integration variable Q' —>■ kQ'. As in [42], we expand the 
integrand around A: = 0, 


( cosh^ — + sinh^ —^ 
\ 2 4 




Cm{x){kQ'' 


2m 


(A.4) 


m=0 


Then the integral over Q' gives the derivative of the delta function: 

r {Q'T = {-2mT6^'^\P - P'). 

7-00 27r 

Thus one can easily perform the integral over P' 

/ QO dp' / P'\^ / P\‘ 

|-(27rA)2™5(2”^)(P-P') ( 2 coshy J = {27ridp)‘^^ hcosh-J 


(A.5) 


(A.6) 


Using these results, we finally get 


Ot^{Q,P) = ^ Cm{x){27rikdp)‘^"^ | 2cosh 


m=0 


P 


= ( 4 cosh^ ^ ~ ^ 


'Kkdp'\j^^^ 


p\P 

2 cosh — 

2 


(A.7) 


B Differential operators 

Here we list the explicit forms of the differential operators up to n = 4. Although 
we have actually computed the differential operators up to n = 17, it is too long to write 
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down and we do not write the explicit forms for n > 5. They are available upon request 
to the authors. 


X){1) 

X){2) 


^2 p2g2a2(i_a2) 

96 (1 +pdn)(l + qd^) ’ 

2 p3g3g3 (l - Qg) ^ ^ ('_1^ _ 12(p+9) ^ 

V 96 y (1 + pdfj,)(3 + pd^)(l + qd^)(3 + 


D(3) = ( — 


^p^g^al+^-gpWip + q)dt 


_ P^Q^dl (1 - (9g) _ 

96 ) (1 +pd^){3 +pd^){5 +pd^){l + qd^)(3 + qd^){5 + qd^) 

+ ^P9 (P^ (89?^ + 464) + 1544pq + 464q'^) ^ (p® i^^q^ + 36) + p^q (93<j^ + 560) 

+ 560pq^ + 36q^'jd^ + ^ (p®q (367g2 + 2272) + 16p^ (7679^ + 336) + 32pq {71q^ + 784) 

+ 5376g^)a)) + ^ (9p^ {23q^ + 8) + p^q {207q^ + 2240) + 448p (Sq^ + 4) + 8g (9g^ + 224)) 9j 
+ ^ {23p^q + 6p2 (3g2 + 7) + pg (23<;2 + 112) + 42g2 + 56) 8^ + ^ (p^ + q^) 


X)W 


(3 ~ ^l) _ 

^ 96/ (1 + pdfj){3 +p9^)(5 + p9^)(7 + p9^)(l + qdf_i){3 + qdfi){5 + g9^i)(7 + 

381p5q5 942 4 4 pV ((3481g2 + 56128) p^ + 147360gp + 561289=) „ 

- 0,1 - p q [p + q)o;, -^ 

1400 175 1400 

- ^P^?^ ((13279= + 4000) p® + 9 (13279= + 22512) p= + 225129=p + 40009®)a® 

- ((143599'^ + 3985289= + 276480) p* + 9609 (11999= + 3776) p® 

+ 649= (62279= + 121204) p= + 36249609®p + 2764809^^) 9^ - ^ ((46039'^ + 183209= + 3200) p® 
+ 9 (46039"* + 1480329= + 106240) p"* + 1929= (7719= + 2749) p® + I69® (11459= + 32988) p= 

+ 1062409'*p + 32009®)^® + (-9 (278599* + 8769929= + 737280) p® 

- 480 (71239* + 320329= + 5760) p* - 649 (137039* + 6752449= + 539136) p® 

- 30729= (50059= + 26036) p= - 1474569® (59= + 234) p - 27648009*) 9® 

- ^ ((54339* + 300809= + 3200) p® + 9 (54339* + 2760489= + 267520) p* 

+ 48 (575I9* + 392769= + 7840) p® + 649 (47O9* + 294579= + 33984) p= 

+ 2569= (10459= + 8496) p + 6409® (59= + 588)) 9* “ (l (209539= + 14880) p® 

+ 30 (4239* + 123769= + 1440) p* + 9 (209539* + 5516809= + 1184256) p® 

+ 48 (77359* + 524769= + 16960) p= + 969 (1559* + 123369= + 20800) p + 9609= (459® + 848)) 9® 

- ^ (5 (9779= + 40) p® + 89 (2019= + 3560) p* + 24 (679* + 19809= + 980) p® 

+ 9 (48859* + 475209= + 139776) p= + 64 (4459* + 21849= + 1200) p 

+ 4O9 (59* + 5889= + 1920)) 9= - (449P® + (99= + 75) p* - 129 (9= - 8) p® 

+ 3 (39* + 609= + 28) p= + (449® + 969® + 4489) p + 759* + 849= + 20o) 9^ - 2304 (p® + 9®) . 


C Derivation of the TBA functional equations for p = 2 

Here we derive the functional equations (2.55). We start with the recursion relation (2.80). 
Defining new functions by {x) = e"^(j)^p {x), then the recursion relation (2.80) is rewrit- 
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ten as 



dy x-y tPe-iiy) 
{2-KkY 2 sinh ^ (2 cosh |)9' 


(C.l) 


Following the argument in [47], these functions also satisfy the following functional relation 


+ 2'Kik) + pp{x - 2-Kik) + 2pp{x) = x-^q PPi{x)- (C.2) 

Therefore the original functions (l)P{x) satisfy 

- (t)P{x + 27Tik) - (pP{x - 2-Kik) + 2(I)P{x) = ^ 2 ^ 0 ^^ x)g 4-\(^)- (C-3) 

Let us introduce a generating functional of pp{x)\ 

00 

$L)(x) = 'PP{—kY pp {x). (C.4) 

1=0 

The functional relation (C.3) is then written as 

+ 2Kik) + — 2Kik) = 2(1 + t{x))^^^\x), (C.5) 

where t{x) is defined by (2.56). We have used the identity: —pp{x + 2Kik) — pp{x — 
2Kik) + 2pp {x) = 0 for PY\x) = x^ {j = 0,1). One notices that the functional relation 
(C.5) is the same form as Baxter’s TQ-relation. The functions <h^-^^(x) (j = 0,1) are two 
independent solutions of the TQ-relation. A crucial fact is that these two solutions satisfy 
the so-called quantum Wronskian relation: 


{x + Kik)^^^^ {x — Kik) — (x — nik)^^^^ {x + Kik) = const. (C.6) 


The constant is fixed by taking the limit «: —)■ 0. Since we have $(^)(x) = 1 -|- 0{k) 
and $(^)(x) = X + 0{k)^ one easily finds that the constant must be —2Kik. For later 
convenience, we rescale 4>*''^^(x) by 

$+(x) = $(°)(x), 4>_(x) =-(C.7) 

Kik 

Then the rescaled functions satisfy the quantum Wronskian 


<l>+(x -|- Kik)^-{x — Kik) — ‘k+(x — Kik)^-{x -|- Kik) = 2. (C.8) 

As shown below, this relation plays a crucial role in deriving (2.55). Our goal is to compute 
the diagonal elements of the resolvent: 

CX) 

^)- (C-9) 

n=0 

Using the formula (2.87), one finds 

R{x) = (4>(0)(x)9,c^«(x) -4>«(x)9,4>(0)(x)) = (C-10) 
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where W[f,g] = f{x)g'{x) — f'{x)g{x) is the standard Wronskian. 

Now we derive (2.55) from the quantum Wronskian (C.8). In the following, we use a 


notation, for simplicity. 


= fix±Trik). 

(C.ll) 

Let us first consider the square of (C.8) 


($+$- _ $-$+)2 = 4. 

(C.12) 


It is easy to see that this is equivalent to 

+ - 1- (C.13) 

Introducing the functions ^{x) and g{x) by 


^(x) = <I>+(x)$_(x), 


1 


i]{x) = -($_|_(x + 7rik)^-{x — nik) + <k+(x — 7riA;)$_(x + 7rik)'j, 

then we get the hrst equation in (2.55). 

Next we rewrite g{x + vrifc) + g{x — irik) as 


(C.14) 


g{x + irik) + g{x — -rrik) = 


1 


<h+(x + 2Trik)^-{x) + <I)+(x)$_(x + 2TTik) 

+ <I>+(x)$_(x — 27Tik) + $+(x — 27riA:)<h_(x) 

= 2(1 +t(x))^(x), 


(C.15) 


where we have used the TQ-relation (C.5). Using r{x) in (2.56), we get the second equation 
in (2.55). 

Finally we consider the Wronskian 


w(x) = -IU[<I>+,$_] = -(<!>+(x)<I>'_(x) — (x)<I>_(x)). 

^ z 

One can see 

w{x + 7:ik) w{x — nik) 1 (4>1|1<I>1)' (<h^$l)' 

^{x + 7rik) ({x — irik) i\_ <I)ljl<I>I 


(C.16) 


(C.17) 


From the quantum Wronskian, we have ($+<I>jl)' = (<I>lj]<!>_)'. Thus one obtains 

w{x + TTik) w{x — nik) 2 

^{x + 7rik) ^{x — irik) i ' 

On the other hand, we have 

g'{x) _ 1 

7/2(x) — 1 ^{x + 7rik)^{x — irik) 2 ' 


(C.18) 


(C.19) 


Therefore we find the final equation in (2.55). It is easy to see the grand potential is written 
as (2.57) by using (C.IO). 
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D Planar solution for p = 1 

In this appendix, we compute the free energy of the (l,g)-model in the’t Hooft limit 


N 

k oo, N —)• oo, A = — = fixed, q = fixed. 

k 

For p = 1, the canonical partition function takes the simple form 

d^x a<,tanh2 2^ 


Z{N,, 


I =^f 

ip=i 2^N\ J 


(27r)^ ni(2cosh^)'? ’ 


(D.l) 


(D.2) 


Zi = , 


which is the one-parameter deformation of the A^j-matrix model by k. If we change the 
variable as 

(D.3) 
(D.4) 
(D.5) 


then we find 


where 


Z{NMr>=i = 


m J ( 27 r)^ rr 




V{z) = qlog{z'^/^ + z-'‘/^). 
In the large-A; limit, this potential is rewritten as 


V {z) = nkq 


^ 1 mLi2(2;)-I-imLi2(-z))-I-0(A: 

Ztt tt^ 


(D. 6 ) 


In [63], the authors have computed the planar free energy of the matrix model with the 
potential 


Vkmz{z) = - 

g 


log z 1 


-k 


27r TT 


— ^ImLi2(i^;e^)-|-ImLi2(z^e ^'^+0{k 


(D.7) 


by using the technique in [64]. Hence if we take g = l/iirkq) and ^ = 7 rz /2 in their planar 
solution, then we can obtain the planar free energy of the ( 1 ,( 7 ) model. Since the ABJM 
case corresponds to g = 1 , this means that the planar free energy of ( 1 ,^) model is the 
same as the one of the ABJM model with the replacement k —)• kq: 


log Z(iV, k) planar = ^Og ^ I (p,q)=l,planar,k^kq ' 


(D. 8 ) 


Recalling that the planar free energy of the ABJM theory has the worldsheet instanton 

_27r ! 

effect with the weight 0{e V *: ) [ 6 , 9], we easily see that the planar free energy of the 
( 1 , q) model has also the non-perturbative effect of the order 


0{e 




(D.9) 


which is the same as the expected WS instanton effect from the gravity side. 
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Let us see that this result is consistent with our result on the grand potential. As 
explained in [2, 26], the ’t Hooft limit in the grand canonical language is 


u —>■ oo, k —)■ oo, A = T = fixed. 

k 

In this limit, we can expand the grand potential as 


(D.IO) 


9=0 


(D.ll) 


which should be considered as the “genus” expansion of the grand potential. Then the 
“planar” grand potential i7o(A) O') is related to the planar free energy Fq{X) by the Legendre 
transformation: 

Jo(A, q) = Fo{X) - A^i"o(A), A = ^^o(A). (D. 12 ) 

Noting the planar free energy takes the form 

Fo(A) = gVo(Ag), with Xq = -, (D.13) 

q 

then the Legendre transformation relation becomes 

= /o(A,) - Xq^FoiXq), ftq = ^ = ^MXg), (D.14) 

This relation leads us to 

Jo(A,g) = 9^^o(A, • (D-15) 

We can easily check that the perturbative grand potential (2.39) and the fist two worldsheet 
instanton corrections (3.23) satisfy this relation: 


and 


lApert (h) ^) — q k 


2 n 1 , 


+ 0{l), 


di(l,g,A:) = + 0(1), d2{l,q,k) 


Sq'^k^ 

IGvr^ 


+ 0 ( 1 ). 


(D.16) 


(D.17) 


E On the g-Gamma function 

In this appendix we propose a useful integral representation of the g-gamma function. The 
g-gamma function rq( 2 :) is defined by 

r,(=) = (1 - (E.i) 

(q ;q)oo 

in terms of the g-Pochhammer symbol 

CX> 

(«;q)oo = 

1=0 
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The two important properties of the g-gamma function are the following functional relation 
and the behavior in the limit q —)■ 1 

rq(z + i) = i^rq(z), iirnr,(2) = r(z). (e.s) 

1 - q ^ q^l 

The infinite product representation (E.l) of the ( 7 -gamma function is well-defined when 
|q| < 1. However, in our case of interest q = with S M, we have to deal with the q- 
gamma function with |q| = 1. In this case, the naive infinite product expression (E.l) per se 
is ill-dehned, and we have to define the ( 7 -gamma function with |q| = 1 as a certain analytic 
continuation from |q| < 1. In the literature, such analytic continuation was proposed by 
using either the double sine function [57] or the Faddeev’s quantum dilogarithm integral 

[65]. 

In this paper, we propose an alternative integral representation of the g-gamma func¬ 
tion with jqj = 1 , which is useful for the numerical calculation of the instanton coefficient 
in (3.15). We regularize the infinite product appearing in ( 7 -Pochhammer symbols by using 
the zeta-function regularization. For q = the ( 7 -Pochhammer symbols in (E.l) can be 
rewritten as 


CO / • h{n-\-z) \ 

(q^;q)oo = n 


CO . h{n+z) 


Sin 


n=0 


n=0 




cx> • h(n-\-l) 


Sin 


h 


n=0 2 

where Cis,a) denotes the Hurwitz zeta function 


n ’ 
2 


(E.4) 


C{s,a) = y- --. 

n {n + ay 
71=0 ^ ' 


(E.5) 


Plugging the value of C{s, a) at s = 0, —1 


C(0,o) = i-a, C(-l,a) = i Q-Fa-a^ 


(E.6) 


into (E.4), we find 


r^(^) = ef(--i)(--2) 


h 

. 2 , 


1-2 


00 „• h{n+l) h 


n 

n=0 


Sin 


h n(n+z) 

o sm —T— 


(E.7) 


Now, let us consider the h expansion of the inhnite product part in (E.7). Using the 
expansion 


log 



2m{2m)\ ’ 


(E.S) 
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we find 


log n 


oo . h(n-\-l) 

sm ^ ^ ^ 


n=0 


sm 


h{n+z) 


oo oo oo 

^ [log(n +1)-log(n + 2 ) +'^Y1 




2m r 


n=0 


n=0 m=l 

t 2 m r 


2m{2m)\ 


(n + 1)"™ - (n + z)"™ (E.9) 


C'(0,1) + C'(0, z) + Y2 ^ ^2m{2m)\ “ C(-2m, z) 


This can be further simplified by using the relation 

r(z) 


C'(0,z) = log 


V^’ 


C,{-2m,z) = - 




2m+l 


(^) 


2m + 1 


(E.IO) 


and (E.9) becomes 


log n 


sm 


h(rL+l) 


n=0 


sm 


h{n+z) 


logr( 2 ;) + 


m=l 


{-irB2mB2m+l{z)h^^' 
2m{2m + 1)! 


(E.ll) 


Putting all together, we find the following representation of the q-gamma function with 

|q| = 1 


1-2 


T^{z) = 64 


_ ^i^( 2 -l)( 2 - 


2 )r(.< 


sm ■ 

W 

. 2 


exp 


E 

_m=l 


{-l)'^B2mB2m+l{z )fi- 
2m{2m + 1)! 


2m 


(E.12) 


Using the property of the Bernoulli polynomial (3.17), one can show that (E.12) indeed 
satisfies the functional relation in (E.3), as required. Also, one can easily see that (E.12) 
reduces to the usual gamma function in the limit —)• 0. 

However, (E.12) is still a formal expression since the summation in the exponential 
factor is a divergent asymptotic series. When 0 < z < 1, we can resum this series by using 
the integral representation of the Bernoulli polynomial (3.18) and (E.8). Einally, we arrive 
at our integral representation of the g-gamma function valid for 0 < z < 1 and |q| = 1 


1 — Z 


rq(2) = eT(^-i)(^-2)r(z) 


sm: 

w 
. 2 


exp 


f 


dt 


sin 2Trz 


cosh 2nt — cos 27r2: 


log 


/ sinh ^ 


(E.13) 


For the case z 0 (0,1), a similar integral representation can be obtained by using the 
functional relation (E.3) repeatedly. 
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F Exact values of Z (iV, k) and instanton corrections 

Using the exact values of the partition functions for various integral (p, q, k), we can deter¬ 
mine the non-perturbative part of the modified grand potential^® 


k) — i/pert(/^) ^)) 


(F.l) 


by the numerical htting, in a similar way as the ABJM case [11]. 

In this appendix, we list the non-perturbative part of the grand potential Jnp(h) k) and 
the exact partition functions Z{N, k) for N = 2,3, 4, for various integral {p, q, k). We drop 
the N = 1 case since we know the exact value of Z{\, k) in a closed form (4.1) for general 
{p,q,k). Actually we have computed the exact partition functions for higher N > 5, but 
they are too lengthy to write down in this appendix. We have also computed the exact 
partition functions for several other {p, q, A:)’s which are not listed below. They are available 
upon request to the authors. 


The case of {p,q) = (1,2). 


-(2,2) ^ .(3,2) . 5". .,4.2) . 


TT^ — 


10247r2 

Jnp(h, 2 ) = + 


vr 


3686407r3 
10p2 + 7p + 7/2 


+ 1 


TT^ 


943718407r4 
-2^, , 88/i-k52/3 3 ^ 

^ 37r 


+ 


269/r2 193^/4-^ 265/16 


vr^ 


-F58 


4792p + 1102/5 

- “T e , 

OTT 


^(2,3) = . ^(3,3) = 


-21384 -k 1331l7r2 - 2048\/37r^ 


Z(4,3) = 


311047r2 ’ ' ’ ' 100 776967 r 3 

614304 - 18213127r2 - 32768V37r3 1962977r4 


193491 76327r4 


8 _2m 


.iM 


"AipChj 3) — e 3 6 e 3 -|- 


4/r^ -|- 2p -|- 1 88 


3tt'^ 




g 2/4 _ 


238 _ 8/4 


9 


-e 3 


848 io<4 
+ 1 ^^' ^ + 


h2p? -h /r -F 9/4 1540 


Z(2,4) = 
Z(4,4) = 


57r2 - 48 


Z(3,4) = 


Ott^ 9 

-2640 -k 8337r2 - 1807r3 


e-^^ + 


82672 _i4/i 
-e 3 , 


189 


81927r2 ’ ^ ^ 58982407r3 

6400 - 157767r2 - 48647r3 308l7r^ 


4026531847r4 


— 2 \/ 2 e 2 -|- 


h + 1 4 

16\/2 3/4 

e"^ H-e"^ -h 

lO/t^ + 7p + 7/2 45' 

TT 

3 

[ 27r2 2 J 


=-2/4 


288^/2 5/4 

+—-—e 2 + 


88p-^52/3 640 

Ott 3 


= -3/4 


^®The modified grand potential J(/4, k) is related to the full grand potential f7(/4, k) by 

= E< 


gj7’(/.t,fc) _ V ^ ^ J(^+27rm,fc) 


As shown in [11], the modified grand potential removes the “oscillatory part” from the full grand potential. 
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Z(2,6) 

Z(4,6) 

JupilJ-i 6 ) 
The case 
Z(2,2) = 

'Aip(/^)2) = 


33l7r^ - 3240 


Z(3,6) = 


-495720 + 2870377r2 - 43520V37r3 


7464967r2 ’ ' ’ ' 24186470407r3 

459794880 - 11613961447r2 - 320716800\/37r3 + 2897742257r‘^ 


8 _M 14 2fi 

of {p,q) = (1,3). 

32 - 37r2 


501530650214407r4 
■2(/i + l) 


37r 


+ 8^3 


154 _4m 

T 


e ^- —e 3 


(F.2) 


122887r2 


7552 - 7657r2 _ 143360 - 2787847r2 + 26 7757r'^ 

“ 17185920 ^’ 


422785843207r4 


2(2u + 3) _2 m 

-6 3 


+ 

+ 

+ 


y/SiT 

{2fi + 3)3 


{2p + 3)^ 


+ 


_4m 
6 3 


76//^ 47(2// + 1) 2(2/i + 3) 


TT'^ 


37/2 


07/2 


V3 


TT 


3-2/1 


_ .§££ 
6 3 


(2//+ 3)^ 2(2//+ 3)2 166//+ 133/4 2 

27/4 ^2 + 3 

\/3(2// + 3)® 8(2//+ 3)3 332//2 + 1129///2 +399/4 4^3(2//+ 3) 

- H--1- n -48 


Z(2,4) = 
Z(4,4) = 

'lnp(/^,4) = 

The case 

Z(2,2) = 

Z(4,2) = 

Jnp{pi2) = 

Z(2,3) = 
Z(4,3) = 
'^np(/^, 3) = 
Z(2,4) = 
Z(4,4) = 


57/3 

1057/2 - 1024 




7/j 


7/^ 


7 / 


10 ^ 

6 3 , 


Z(3,4) = 


-1024 - 26247/+ 9397/2 


7864327/2 ’ ^ ^ 1006632967/2 

367001600 - 8713953287/2 - 3513753607/3 + 1963694257/4 
4329327034368007/4 ^ 


+ 

1 

[2(2/7+ 3) , 1 

2m 

e 3 + 

3V3(2/7 + 3) 


V37/ 


TT 




(F.3) 


of (p,q) = (1,4) 

277/2 - 256 


Z(3,2) = 


219757/2 - 216832 
3932167/2’ ^ -j - 130799370247/3 ’ 

3153920000 - 70923371527/2 + 6862259257/4 
160726266150912007/4 ^ 


\/2(2// + 4) _M 
- — - -e 2 + 


TT 


(2// + 4)2 3(// + 1) 

O “T 


TT^ 


TT 




+ 


2^2(2//+ 4)3 
37/3 


-4\/2 


_3^ 

6 2 , 


1791557/2 - 1767096 


Z(3,3) = 


-585293688 + 9080285497/2 - 155975680^37/3 


1813985287/2 ’ ' ’ ' 74054135070727/3 

1168603329600 - 27060642557287/2 + 120082923520^37/3 + 1959795867757/4 


16 _E 86 


284367878671564807/4 
-^^ + 88 ^ 3 ^-" 2458 .4^ 


37 / 


e ^- T^—e 3 , 


58417/2 _ 57344 


Z(3,4) = 


-14508032 + 56827117/2 - 13409557/3 


1509949447/2 ’ ' ’ ' 33484638781447/3 

1192021196800 - 25696642744327/2 - 8199206830087/3 + 5091130166857/4 

167576127384689049607/4 
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T / A\ ^ —ii 

<Aip(/^)4) = —: —e 4 -|- 


sin ■ 


V2(// + 2) 


+ const 


— E 

e 2 . 


Trtan ■ 

The case of {p,q) = (1,6). 

n^ 33757r2 - 32768 1372161757r2 - 1354203136 

Z(2,2) — 9 —, Z(3,2) — 


(F.4) 


Z(4,2) = 


7549747207r2 ’ 66072367595520007r3 

18673845640626176 - 364554470129664007r2 + 35020036770468757r4 


T / o^ _ 2 ( 2 /i + 6 ) _E 


203688782836089539788800007r4 
(2// + 6 )" ^ V3(4/i + 6 ) _ 2 


TT 


TT^ 


TT 


2/j, 

e ^ 


+ 


( 2 //+ 6)3 3^3(2//+ 6 ) (4//+ 6 ) ^ 32//+ 56 8 




27r2 


37r 


V3 




(F.5) 


The case of {p,q) = (2,3). 


Z(2 1) = ^ - 48) z,3 n _ 6784 - 687^^ 

^ ^ 409607r2 ’ ^ ^ 707788807r3 ’ 

45731840 - 923519287 r 2 + 88877257r4 


Z(4,l) = 


T f n _ 4(2/^ + 3) _ 2 m 


TT 


104639496192007r4 

(2/i + 3)2 3v^(4/i + 3) 10 

_ 


vr^ 


27r 


+ 


y/STT^ 


+ 


2 ( 2 //+ 3)3 ^ 8(4//2 + 49///2 + 49/4) , 4(2/i + 3) 
^^^2 + VSvr 

-240 + 18707r2 - 1877r'^ 


_4ii 

e 3 


+ 16 


= -2/4 


, , 457r2 — 436 , , 

0 . 044444,2 > ^(3,2) = 


Z(4,2) = 


9830407r2 ’ ^ ^ 2359296007/5 
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234392471470080007/8 

«^np(l^,2) = 4\/3e 3 + 

3//2 + 9// - 109/3 67// + 17 391 


^/3{2^^ + 3) 37 

2 /j, 
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r 8// + 20^22V3l 

TT 3 


TT 
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+ 


_4ii 

e 3 , 


"s/St/ ' 

-6076 - 21607/ + 13057/2 
314572807/2 ’ 

-15360 - 576007/ + 565607/2 9562807/3 + 4319727/“^ - 2355757/5 


Z(2,4) = 

“ 4831838208007/5 

Z(4,4) = (7096320000 - 336211814407/2 -495500544007/3 +238227230780//'^ 

+14424390288007/3 - 18920352743527/® - 20499105904207/^ 
+7954396931257/8)/(240017890785361920007/8), 

Jnp(//,4) = 12 ^/ 2 e-t + (-26 + 5^/3)e-t. 
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The case of {p,q) = (2,4). 


Z{2,1) 
Z(4,l) 
Jnpip: 1) 
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Z(4,2) 
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Z(3,3) 
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-10395 + 99007r2 - 8967r4 
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TT 


3347646382080007r8 
2^2 _ 18^ - 18 
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64(/2 + 2)^ -295/i2/3 + 3767/i/18 - 4945/36 
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-9461760 - 11623040772 + 55137032774 - 545737577® 
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(-6125543424000 + 15615912632320772 + 712581545728077"^ 

-1372544905622477® + 130192775212577® )/(70205233054718361600077® 

8V2e-U -22 e-U - ^ ^ 


22 ] 

6 2 + 

16^2(3/7 + 4) 373 

TT 


ro 

_ 1 


6200145 - 5356836772 _ 752540 ^ 377 ® + 896000774 
85710804480774 ’ 

(237546155385 - 79588879722077^ - 323485747200^377® 

+175886392204877"^ + 353173708800^377® - 34693120000077®)/(222690040924538880077®), 
(-60294434727802275 + 23227323309194004077^ - 133698981198182400^377® 
+49864332743345016077"^ + 393619532866560000^377® - 102396466118698444877® 

-238080482537963520^377^ + 20842428006400000077® )/(9454742159931510246604800077®), 

32 . , 

—pC 6 — 16e 3 . (F.7) 

v3 


The case of {p,q) = (2,6). 


Z(2,l) 

Z(4,l) 


JnpiPi 1) 


Z(2,2) 

Z(3,2) 

Z(4,2) 


93772 _ 77 Q -18293275 + 12655344772 - 1094400774 
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-39855355314775 + 64469243716878772 - 24750110932272774 + 188732643840077® 



14772580993648558080000778 

8(/7 + 3 ) 

4(/7 + 3)2 3V3(2/7 + 3) 10 

2pL 

-F-e ® + 

■s/377 

772 ^ 77 3 

e 3 


, [16(/i + 3)® 8(4/i2 + 49/i + 49) , 8(/i + 3) 

+ ---^-F-Id ® I 

■s/377® 3772 ■s/377 
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20927899238400774 ’ 

-51211183063040 + 7706727614976772 + 52775122747500774 - 537308584312577® 

255747634699331174400077® 

(-25823407803062655385600 + 50663329067554514927616772 
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+186487951690759002877447r^ - 396457161224434761520807r® 
+37755261796563792731257r®)/(5075820179618539253692524134400007r®), 
10(/i + 3) 


'Aip(/^)2) — 24e 6 -|- 

vSvr 

The case of {p,q) = (3,4). 
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_i£ 

e 3. 


(F.8) 


Z(2,l) 

Z(3,l) 

Z(4,l) 


JnpifJ-i 1 ) 


Z(2,2) 

Z(3,2) 

Z(4,2) 


2) 


-2880 + 136807r2 - 133367r‘^ + 12157r6 
88473607r6 ’ 

241920 + 10281607r2 - 36318247r'^ + 25548327r® - 2226697r8 
10701 76665607r9 ’ 

(-7770470400 - 270768960007r2 - 5674252953607r^ + 13040549508807r® 

-7647020842327r® + 646862658757r^°)/(8438128972922880007r^°), 


2^3(2/!+ 3) 915 

2/j, 

e"” + 

r 28/i + 34^^g^^1 

71 18 


TT 


-241920 + 7812007r2 - 13318967r^ + 1275757r® 

39636 1 728007r6 ’ 

-10866240 - 662244007r2 + 1029422687r^ - 97390657r® 
6696927756288007r7 ’ 

(16428264652800 - 2654804984832007r2 + 9708669479308807r‘^ 

-321512581 74963207r® + 56967892184141127r® - 35472971495340487r^° 


+3041779559568757r^2)/(20967962939009217331200007r 


12 ' 


8\/6e 
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(F.9) 


The case of {p,q) = (4,4). 


Z(2,l) 

Z(3,l) 

Z(4,l) 


JnpifJ-i 1 ) 


-7875 + 396907r2 - 2327571^ + 203271® 

406425600718 ’ 

675675 - 11434507i2 + 1703740571^ - 25454 77071® + 995491271® - 76492871^® 

212446789632007112 ’ 

(82084377375 + 1057890023190071^ - 907551930285071^ + 8429530663554071® 
-26455473328784571® + 22682227421089671^® - 6533438552409671^^ 


+45575850147847i^'^)/(304002178262079897600007i1®), 


12/1 + 28 

6 2 + 

64(/i + 3) , 736' 

TT 


71 3 


+ 


19/i2 + 127/1 512/1 
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+ const 


TT^ 


TT 


0-1^ 


3/i. 

4 


(F.IO) 


The case of {p,q) = (4,6). 

^ -872025 + 21101857i2 - 103488071^ + 8492871® 

^ ^ “ 178827264000718 ’ 

Z(3,1) = (1355404050 - 185966030257i2 + 8499448503071“^ - 8771899292071® 

+2833634419271® - 20570112007i^®)/(30731393787494400007i^2)^ 
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Z(4,1) = (-25947408325460997375 + 980597970505170517507r2 - 1595964317417656436257r‘^ 
+13673391282168341724007r® - 27442524814765174514407r® 
+18780376465788954016007r^° - 4764961742101140602887ri2 
+3 1 7112 1 74820196352007r^'^)/(476980876903659017317908480000007r^®), 


Jnp(//, 1) = 16\/3e 6 + 


4(/i + 3) 
^/3^T 


- 50 


_i£ 

e 3. 


(F.ll) 


The case of {p,q,k) = (l,q, 2). In the case of {p,q,k) = {l,q,2) with general q, from 
the numerical fitting we find that the non-perturbative part of grand potential takes the 
similar form as the A^j-matrix model [28] 


JnpiPj k — 2 ) — ^ ^ m{,P: 

m=l 


q 


+E 

i=i 


+ ci{q) 




+ 


(F.12) 


where Vm{p,q) is a order polynomial of p and the ellipses denote the contributions of 
bound states. The first three terms of Vm{p, q) are given by 


'Pi{p,q) = —Pu 

Si 


P2{p,q) = -Pi + 


3S3 


Pz{p,q) = 
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3si 


2p3_ 
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where we defined 


vrn 2np, + q 

Sn = sm —, Fn = - 

q TT 


(F.13) 


(F.14) 


We also conjecture that the 1-instanton coefficients bi{q),ci{q) in (F.12) are given by 


bi{q) = 2ci{q) = 


1 r^(-g) 

271^ cos(7r(7)F(—2(7) 


(F.15) 


By taking the limit q ^ n (n £ Z), one can check that the conjectured form of instanton 
coefficients (F.13) and (F.15) correctly reproduces the result of Jnp(//, 2) listed above for 
the {p,q,k) = {l,q,2) case with various integer q. 

One can derive the expression of bi{q) in (F.15) by taking the limit A: —)■ 2 of 'yi{q, k) 
given by (3.9). However, if we take the limit naively, we get a wrong answer. To reproduce 
(F.15), we have to first rewrite 7i(g,A:) by using the transformation of hyper geometric 
function as 


ii{q,k) = 


1 r2(-g) p 

7r2A:F(—2g)^ ^ 

1 r^(-g) . 

tt^A; F(—2(7) ^ ^ 


( <1 . 

q 

1 

V 2 ’ 

2’ 

2 


2 '^k\ 
g;sm -J 

1 o'^k 


(F.16) 


Taking the limit A: —)• 2 in the last expression, we correctly obtain bi{q) in (F.15). 
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